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NOTATION 


1.  Unit  means  service  unit  or  counter  at  which  the  server 
is  serving. 

2.  L.S.T.  is  abbreviated  for  Laplace-Steiltjes  Transform. 

3.  H(x)  is  the  distribution  function  of  service  times. 

Its  L.S.T.  is  denoted  by  h(s)  and  its  first  three 
moments  by  a,  0  and  Y  respectively. 

4.  Y.(s)  denotes  the  L.S.T.  of  the  distribution  of 
busy  period  of  an  m|g|1  queue  with  service  time 
distribution  H^(x). 

5.  Y(s)  denotes  the  L.S.T.  of  the  distribution  of  busy 
period  of  the  whole  system  (All  the  service  units 
considered  together). 

6.  The  convolution  of  two  distribution  functions  F(x) 
and  G(x),  0  <  x  <  is  denoted  by: 

F*G(x)  =  [  F(x-u)  dG(u) 
o 

The  m-fold  convolution  of  F  is  denoted  by: 

F^(x)  =  F*F^m’1^(x) 

7.  U(x)  is  the  unit  distribution: 

U(x)  -  0  if  x  <  0 
=  1  i  f  x  >  0 

8.  6..,  is  the  Kronecker  delta  defined  by: 

id 

5  .  «  0  if  i  /  j 
*  1  if  i  -  J 


For  referring  the  equations  we  use  the  following 
convention:  (n)  means  the  n-th  equation  of  the 
present  chapter  and  (m*n)  means  the  n-th  equation 
of  the  m-th  chapter. 


9. 
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CHAPTER  I 

A  SINGLE  SERVER  TANDEM  QUEUE  WITH 
NON-ZERO  SWITCHING  IN  UNIT  1 


1.  Concepts  and  Definition* 

In  this  chapter  ve  consider  a  queueing  process  with  two 
service  units,  unit  1  and  unit  2,  and  a  single  server.  The 
server  attends  to  the  two  units  alternately  according  to  some 
switching  rule.  A  switching  rule  [Neuts  and  Yadin,  1968]  is  a 
rule  describing  how  the  server  changes  from  one  unit  to  the 
other.  The  server  may  change  from  one  unit  to  the  other  either 
by  a  non-zero  switching  rule  or  by  a  zero  switching  rule.  By 
a  non-zero  switching  rule  the  server  continues  to  serve  in  a 
unit  until  some  upper  number  of  consecutive  services  have  been 
cocipleted  and  then  he  switches  to  the  other  unit.  By  a  zero 
switching  rule  the  server  stays  in  a  unit  until  the  queue  in 
it  becomes  empty  and  then  he  switches  to  the  other  unit. 

In  this  chapter  we  discuss  a  non-zero  switching  rule  for 
unit  1  and  zero  switching  rule  for  unit  2.  The  zero  switching 
rule  for  unit  1  is  dealt  in  the  next  chapter. 

We  say  that  two  units  are  in  tandem  ^h %n  the  output  of  the 
first  unit  is  the  input  to  tne  second.  It  is  assumed  that 
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customers  arrive  in  unit  1  in  accordance  with  a  Poisson  process 
of  density  X.  The  input  for  unit  2  is  those  who  have  completed 
service  in  unit  1. 

The  durations  of  the  successive  service  times  in  units  1 
and  2  are  identically  distributed  independent  positive  random 
variables  with  distribution  functions  H^(*)  and  Hg(‘) 
respectively.  Further  the  service  times  are  independent  of  the 
arrival  times. 

In  the  case  of  non-zero  switching  the  server  starts  in 
unit  1  at  time  t  =  0  and  continues  to  serve  in  it  until  he  has 
given  k  services  without  interruption  or  until  the  queue  be¬ 
comes  empty,  whichever,  comes  first,  k  is  a  positive  integer  which 
we  will  call  .  as  the  switching  parameter.  The  time  interval 
epent  without  interruption  in  unit  1  is  called  a  1-task. 

Similarly  we  define  a  2-task.  A  1-task  followed  by  a  2-task 
both  together  will  be  called  as  a  cycle  of  tasks. 

The  customers  who  have  completed  service  in  unit  1  queue 
up  in  ftront  of  unit  2.  The  e«rver  after  completing  the  1- 
task  switches  to  unit  2  and  nerves  there  until  the  queue  in  it 
becomes  empty.  After  finish tr>g  the  task  in  unit  2  the  server 
switches  back  to  unit  \  and  continues  the  process. 

When  k«l,  we  obtain  sl&pir  an  M/o/l  queue  with  service 
time  distribution  *  Hg( • } 
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2.  Distribution  of  Busy  Period 
The  server  bejins  in  unit  1  at  t«0  and  serves  between  the 
two  units  alternately  according  to  some  switching  rule.  The 
tine  required  for  both  the  units  to  become  empty  simultaneously 
for  the  first  tine  is  called  a  busy  period  of  the  systen. 

Suppose  that  there  is  a  Poisson  input  of  density  X  in 
unit  1  and  that  the  service  time  distributions  of  the  two  units 
are  H^C*)  and  Hg(‘)‘  Then  since  the  distribution  of  busy 
period  does  not  depend  upon  the  order  in  which  the  customers 
are  served  [Welch  (1965) ]>  the  distribution  of  busy  period  of 
the  model  defined  above  is  equivalent  to  the  distribution  of 
busy  period  of  an  M/G/l  queue  with  input  rate  X  and  service 
time  distribution  the  convolution  *  Hg( • )•  Hence  from  the 
classical  results  of  an  M/G/l  queue  [Takacs  1962,  p.  4 f]  that 
if  y(s)  is  the  Laplace  Stieltjes  Transform  (L.8.T.)  of  the 
distribution  of  busy  period  then  y(s)  is  the  unique  root 
in  the  unit  disk  jz|  <  1  of  the  equation 

(1)  z  »  hj^s  +  X  -  Xs)  hg(s  +  X  -  Xz), 

The  expected  length  of  busy  period  is  given  by: 

(2>  -  V'(0+)  ■  if  1-Xo^-Xo.jX), 


if  l-x  aj.  -  \  Oj  -  o, 
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3.  The  Basic  Imbedded  Semi-Markov  Process 
and  its  Transition  Probabilities 
Macros copically  the  queueing  process  consists  of  busy 
periods  alternating  with  idle  periods.  Each  busy  period 
consists  of  a  random  number  of  alternating  1-tasks  and 
2-tasks.  Every  busy  period  can  be  decomposed  into  a  random 
number  of  cycles  of  tasks. 

Hare  we  assume  that  at  t=0  there  are  i  >  0  customers 
in  unit  1  and  none  in  unit  2.  In  the  case  i=0  the  process 
starts  with  an  idle  period. 

Let  us  define  the  sequence  of  random  variables  TQ,  T^, 

T9> . .  where  T  =  0  and  T  is  the  duration  of  the  nth  cycle 

c  on 

of  tasks,  n=l,2,....  Let  ^  denote  the  number  of  customers 
in  the  system  at  the  end  of  the  nth  cycle,  n=l,2,....  and 

So’1- 

It  follows  from  the  regenerative  properties  of  the  input 
and  service  processes  that  the  bivariate  sequence  of  random 
variables 

(3)  (v  V  "  2  °} 

is  a  Semi-Markov  sequence  with  state  space  (0,1,2, . }. 

We  recall  the  definition  of  Semi-Markov  sequence. 

Consider  a  double  sequence  of  random  variables 

defined  on  a  complete  probability 

space  and  such  that: 
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U)  W 

p{xo  -  0]  -  1, 

(ii) 

p(jo-k]«a^,  where  a] 

(iii) 

PlV“.  X„  i  x  1 

apace 


=  P{Jn»k,  X <  x  |  J  )  *  Q  (x), 

n-1 

for  n*l,2, ..... 

then  the  process  t(Jn»Xn)>  n  >  0)  is  called  a  Seal -Markov 

sequence.  The  functions  (^(x),  i>j*l,2, _ are  mass  functions 

which  are  non-decreasing  and  they  satisfy: 

Q^j(x)  =  0  for  x  <  0, 

OyH  =  i  J=1.2. . . . 

where  (P^)  is  the  transition  matrix  of  the  Markov  chain 

tJn>  n  >  0}.  For  iturther  details  of  Semi-Markov  sequences 

we  refer  to  Pyke  (1961),  Neuts  (1966). 

To  study  the  transition  probabilities  of  the  Semi-Markov 

sequence  we  first  define  an  auxiliary  probability  function 
« (h)/  \ 

°u  >x)- 

Let  us  define: 

(5a)  g[°'(x)  31  8ij  °(x)* 

where  &  is  the  Kronecker  delta  and  U(  )  is  the  distribution 
Regenerate  at  zero.  For  n  >  1,  G^'(x)  is  the  probability 
that,  in  an  m)g]1  queue  of  input  rate  X  and  service  time 
distribution  H^*)  the  initial  busy  period  involves  at 
least  n  services,  that  the  n-th  service  is  cosrpleted  before 


time  x  and  that  at  the  end  of  the  n-th  service  there  are  j 


customers  waiting,  given  that  there  were  i  customers  initially. 
Then  for  i  >  1: 


(5b)  gW(.)  =  f  e'Xy  ^(y)  , 


(5c>  =  't  f  oWfx-yjr^  , 


v*l  0 


Let  gf^(s)  he  the  L.S.T.  of  G^(x)  and: 
i  j 


n  >  1, 


sjn^(s.z)  “  S  g^(s)zj  ,  |*J  <  1,  n  >  0  , 


Then: 


g[0)(s,z)  =  z1  , 

g^(s,z)  =  z1'1  h^s+X-Xs)  , 


g(nn;(t>t) ,  >[gu 


[g[n)(s,z)  -  g^(s,0)^j,  n>l, 


where  h1(*)  ia  the  L.S.T.  of 
Successive  substitution  yields: 
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(8b) 


•  *  1  , 

From  the  definition  of  G^(x)  it  follows  that: 

(9)  g|n^(s,0)  =0  for  i  >  n  , 

Hence  from  (8): 

(10)  g[n^(s,z)  =  si_n  h”(s+\-\2)  for  i  >  n  >  0  , 

A  Summary  of  Known  Results 

(n) 

The  properties  of  the  probability  functions  (x), 
already  known,  may  be  summarized  as  follows:  For  proofs 
of  these  we  refer  to  Takacs  (l960),Neut3  (1968b) 

Lemma  1.1 

If  G1(x)  is  the  distribution  of  busy  periods  for  an 
M|G|i  queue  with  input  rate  X  and  service  time  distribution 
H^( • )  and  Y^(s)  its  L.S.T.,  then: 

(11)  2  B^i)  *  Y*(s)  >  i  >  1, 
n=l 
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Lemma  1.2 

If  Y,(«,u)  -  I  sjn)(s,0)  c"  ,  |o|  <1  , 
n=l  i 

then: 

(12)  E  4n)(s>0)  on  =  vJ(»,c)  ,  1  >1  , 

n=l  A 

For  ia=  1,  y^(s,1)  =  Y^(s)  and  then  lemma  1.2  reduces 
to  lemma  1.1. 

Lemma  1.3 

If  |z|  <  1,  |uj  <  1  and  i  >  1  then: 


(13) 


„  (n),  .  j  n 

*  £  g;/(s)zd  c 

n=0 


z[zi-V^(s,c)] 

z-ah^a+X-Xz) 


) 


For  o=l  one  may  rewrite  (13)  as: 


(1*0 


00  CO 
£  £ 
n=0  j=l 


z  -  h^s+X-Xz) 


J 


Lemma  1.4 

If  R(s)  >0  and  jo)  <1  then  z  =  y^(s>w)  is  a  root  of 
Takac**  functional  equation: 

(15)  z  =  o  h1(a+X-Xz)  ,  |  z  |  <  1  , 

Fxirther  z  =  y^(B}d)  is  the  only  root  of  this  equation  in  the 

unit  circle  |z  |  <  1  if  R(s)  >  0  and  |oj  <  1  or  R(s)  >  0  and 

M  <1  or  R(s)  >0,  Ju|  <  1  and  1  -  X  ^  <  0. 
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Lemma  1.5 

In  lemma  1.4  tailing  u  =  1  we  get  that  for  R(s)  >  0, 
z  *  Y^Cs)  is  a  root  of  the  equation: 

(16)  z  =  h^s+X-Xz)  , 

Further  8  *  Y^(0)  is  the  smallest  positive  real  root  of  the 
equation: 

(17)  0  =  hL(X  -  xe)  , 

and  if  1-X  <  0  then  9  <  1  and  if  1  -  X  >  0  then  0=1. 

From  lemmas  1.4  and  1.5  it  follows  that: 

(18)  Y^(0  +)  =  3^^  if  i-x«l>o  , 

=  •  if  1  -  X  03  =  o, 

,,  01 

(19)  Y,  (0+)  =  - i - *  if  1  -  X  cr.  >  0  , 

1  (1-X  a,)3  ^ 

If  Y-^fOjo)  *  f(o)  then 

(20)  f '(1)  *  3—-^-  if  1  -  X  03  >  0  , 

=  ®  if  1  -  X  0^  *  0, 
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2  X  a. 


*2  3, 


(21)  f'(l)  - - ^  if  1-X  a,  >  o, 

(l-xc^r  d-x^)3  1 

Now  we  define  the  transition  probabilities  of  the  Semi- 
Markov  sequence  Tn,  n  >  0]  defined  in  (3)  as: 

(22)  \,W  -  P(^  =  3,  I„<x  1^-  i> 

Par  i  >  0  and  j  >  0, 

(23) 

sjw  -  Kf  £  <'<“^'x(v'u)  ^ J'v  ^  4k,<~> 


v=0  o  "u 
(u)(v) 

k-1 


£  £  f  ao<;)!«>)e-X(v'u)  JMv^2lI_  ^  „(r)(v.u)i 

(*)(v) 

where  H^(>)  is  the  n-fold  convolution  of  H(*)>  The  second 
term  on  the  right  hand  side  of  the  above  expression  vanishes 
for  i  >  k. 

If  q^s)  is  the  L.S.T.  of  Q^(x)  and 

*  1 

(24)  q1(s,z)  =  £  q  ,(s)  z'  ,  |*!<1, 

j=0 

then: 

(25) 

,  *  1  (k),  ,  r”  -fm)>  (U)3"v  ,  „(k),  , 

^  ‘  Jo  g‘v (,)  J0  e  Ij^yr  11  h2  (x) 


"l1  gW(.)  f  d  HW(x) , 

r*i  10  Jo  d 


for  i  >  0  and  j  >  0, 


! 

i 

) 
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(26a) 


^(sjz)  =  g^(s,z)hg(s+X-\z)  +  £  g^r^(s,0)hg(s+\-Xz), 


k-1 


k  S  J 

for  i  <  k, 


(26b)  =  g^(s,z)  hg(s+X-Xz),  for  1  >  k  , 

From  (9)  and  (26)  it  follows  that.; 

(27)  qi(s,0)  =  Z  g[r)(fl,0)^(s+X)  if  i  <k, 

A  r=i 


*  0 


if  i  >  k  . 


Next  we  introduce  the  taboo  probabilities^  '(,<) 


defined  by: 

(28) 

o^W  =  P{To+Tl+‘"+Tn  ix>Sn=j>?\/0  for 

I  ^0  =  i}  ’  n  >  1  » 

oifw  -  6ij  ^ 

That  is,  oQ^(x)  is  the  probability  that  a  busy  period  has 
at  least  n  cycles  of  tasks,  that  the  n-th  cycle  ends  not  later 
than  x  and  when  it  does  j  customers  are  waiting,  given  that  i 
customers  were  in  unit  1  at  t*0. 

From  the  definition  it  follows  that: 


(29) 
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and 


■  3\  r  «  «vj(u)- n  2 1  • 

v=l  0 


Since  each  cycle  of  tasks  can  have  atrnost  it  services  in  units 
1  and  2,  in  the  above  formula  v  can  be  atmost  j+k. 


Lrt  A? 

(s)  be  the  L.S.T.  of  qQ^I 

(30a) 

or«(,)  •  io  ^”1(s)- 

(30b) 

o\j(s>  ■  j,  o’«(s)  > 

(31a) 

ori(a>s)  =  J?0  orij(a)  z°’ 

(31b) 

oV*’*)  =  ^ 

(32) 

^C.')  ■  j0 

Then: 

or  |z  |  <  1,  R(s)  >  0, 


length  of  time  a  (virtual)  customer  arriving  at  t  has  to  wait 
before  beginning  service  in  unit  1.  For  the  non-sero  switching 
case  the  virtual  waiting  time  at  time  t  will  be  denoted  by 
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(k) 

'(t),  where  k  is  the  switching  parameter  defined  earlier. 
Queue  length: 

The  number  of  individuals  in  the  system  at  time  t  still 
requiring  some  service  in  unit  1  is  defined  as  the  queue  length 
at  time  t.  For  the  non-zero  switching  case  this  quantity  will 
be  denoted  by  ^k^(t). 

Let  6  .(t,x)  be  the  joint  distribution  of  the  queue  length 
5^(t)  and  virtual  waiting  time  ^k^(t),  given  that  at  t=Q 
there  are  i  >  1  customers.  That  is: 

(36)  e^t.x)  3  *  j,  ^\t)  <  x  |  5(k)(0)  »  i], 

Further  let  for  i  >  1, 

(37) 

$ij(t,x)  =  p[^k\t)=j,  0  <T^(t)  <*,  T|^(t)^0 
for  all  x  €  (0,t]  j  §^(0)  =  i]  , 

Formula  (36)  can  be  written  in  terms  of  (37)  as: 

(38) 

rt 

^(t.x)  -  ^(t.x)  r  J  $.^(t  -  u,  x)  d  M^(u) 

+  p{^k)(tH,  n[k)(t)=c  I  c(k)(0)  *  i]  U(x)  , 


15 


where  M^(*)  is  the  renewal  function  of  the  general  renewal 

process  formed  by  the  beginnings  of  busy  periods, 

.  0  for  x  <  0 
U(x)  *  { 

1  for  x  >  0 

To  obtain  the  equation  (38),  consider  the  event  on  the  right 
hand  side  of  (36)  which  can  be  split  into  three  mutually 
exclusive  events: 

(i)  The  time  t  falls  in  the  initial  busy  period, 

5^(t)  =  j  and  0  <  ^k)(t)  <  x,  given  that  ^k^(0)  =  i. 

(ii)  The  time  t  does  not  fall  in  the  initial  busy  period  but 

in  some  other  busy  period  which  started  at  time  u  (0  <u  <t) 

with  a  single  customer,  $^(t)  =  j  and  0  <  T|^(t)  <  x, 

(k) 

given  that  £  (-)  =  i. 

(iii)  The  server  is  idle  at  time  t  (that  is  H^(t)  *  0), 
given  that  5^(0)  =  i. 

The  probabilities  of  these  three  events  give  respectively 
the  three  terms  on  the  right  hand  side  of  (38). 

For  i  >  1,  let  .  (t,x)  be  the  probability  that  at  t 
the  original  cycle  of  task  has  not  yet  ended  and  that 
^(t)  =  j,  0  <  Tl[k^(t)  <x  and  T^(t)j<0  for  all  T€(0,i], 
giv'..  that  at  t=0  the  service  started  in  unit  1  with  i 
customers.  Then: 

(39)  Vj(t,x)  -Elf  d0^"}(u)  )  (t-u,x)  , 

J  n=0  v=l  o  J 

(u) 
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This  formula  is  obtained  from  the  fact  that  at  time  t  the 
server  is  serving  in  the  (n+l)th  cycle  (n=0,l,,...)  of  the 
initial  busy  period.  The  n-th  cycle  ended  between  u  and 
u+du  (0  <  u  <  t)  leaving  v(=l,2,...)  customers  in  the  system. 

We  define  the  following  transforms  for  R(s)  >  0, 

R(0  >  0  and  jz  |  <  Is 

CD 

0*j(t,s)  =  J  e‘8x  d  e^Ct,*)  > 

«  \  e-^'  0*,(t,s)  dt, 

00 

-i 

=  Z  9  ( Z, a)  zJ  , 

j=0 

■  I0  e'sx  a  ’ 

**.  .  r“  -Ft  *  /  . 

*  Jo  e  Cij(t’s)  dt  ’ 

**(*»■»*)  ^  Z  §**(<,, s)zJ  , 

1  j=0  1J 

^ij(t,a)  =  r  e"SX  d  ',;ij^t,X^  * 

***($>«)  ■  J"  e"?t  dt> 

O 

1  0*0 

<D 

B»l(s)  *  J0  d  K^t)  , 
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Lemma  1.6 

For  R(s)  >0,  R(^)  >  0  and  {2  j  <  1,  the  transform 

is  given  by: 

(40) 

k„i 

£  E  [hg  1(s)C^h2(s)-l]D+x-s-\um*h1(8)h2(s)] 

m=0 

-i-l 

.  [  ah1{e)-h1(fc>+X-\omzh1(B)h2(s))]  J 

•  [<’i1^(s)-l}[z[h1(s)-h1(5+X-Xonzh1(s)h2(s))][anjzh1(s)h2(B)] 

-  [zh1(s)-h1(i+X-Xti)mzh1(8)h2(s))l  qi(|1o#  zh1(8)h2(s)) 

+  (z-l)h1(b+x-Xc.mzh1(s)h2(6))h2(s)g|1C/(tl,tym.;h1(s)h2(8)) 

+  Z  C*(l-c^)h1(s)+(c^z-l)h;L(t>+X-Xeinzhi(8)h2(8))3h2(8)^(i>0)]l  , 
v=l  * 

where  c:Q,  . .  .  .  are  the  k-th  roots  of  unity. 

Proof: 

The  probability  *^j(t,x)  is  given  in  terms  of  the 
probabilities  G^(u)  by; 
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(H) 


J  nt 


£J  i 

V=0  0 

(u)(v)  (vx) 


•t+xrt+5Wk) 

t  Jv 


iv 


(u) 


.-X(t-u) 

(FvTT 


(y  \  (4\  ([y]k) 

dvH^;(v-u)dv^(H^^H2  *  )  (vL-v) 


"i1  f  fT*  ®<r)(u)  Efesui 

r.i  Jc  Jt  J»  10  ” 

(“)(')  (vx) 


dvH^r  )(v-u)d^(H^^H2'  )(vrv) 


,(j), 


( tjw 


+  k-1  J  t  t+x  t+x  ,p)  x(t^j 

r=0  v=l  o  Jv  iv  (FvTT 

(u)(v)  (V]) 


J-v 


(rAi£]k) 

•  d^Cv-uJd  k  )  ,(vv)  * 

where  [  ^  ]  is  the  greatest  integer  not  exceeding  ^  . 

The_firet_term  is  obtained  by  assuming  that  the  server 
is  performing  a  2-task  at  t  and  that  the  previous  1-task 
consisted  of  k  services.  The  cycle  of  tasks  in  which  the 
server  is  serving  at  t  started  with  i  customers  in  unit  1 
at  t«0.  The  1-task  ends  after  k  services  leaving  v 
customers  in  unit  1  at  time  u.  The  number  of  arrivals 
between  times  u  and  t  is  j-v  so  that  at  t  there  are  j 


I 
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customers  in  unit  1.  The  services  of  the  j  customers  start 
after  the  completion  of  the  2- task  in  progress  at  t.  Let 
this  2- task  end  at  time  v.  If  j  is  a  multiple  of  k,  say  mk, 
where  m  is  a  positive  integer,  then  the  (virtual)  customer 
eneters  service  after  the  service  completion  of  j  customers 
in  both  the  units.  That  is,  after  the  completion  of  m  cycles 
of  tasks,  each  cycle  consisting  of  k  services  in  each  unit. 

If  j  is  not  a  multiple  of  k,  say  mk+r  (0  <  r  <  k),  then  the 
(virtual)  customer  enters  service  after  the  completion  of  m 
cycles  of  tasks  together  with  a  further  service  completion 
of  r  customers  in  unit  1.  Let  this  service  completion  occur 
at  time  V^.  Now  we  integrate  and  sum  over  all  choices  of 
v,  u,  v,  Vj_. 

The  second  term  is  obtained  by  assuming  again  that  the 
server  is  performing  a  2-task  at  t  and  that  the  previous 
1-task  consisted  of  only  r  (  <  k)  services,  leaving  the 
unit  1  empty  at  the  end  of  the  1-task  at  time  u. 

The  last  term  is  obtained  by  assuming  that  the  server 
is  performing  a  1-task  at  t.  The  cycle  of  tasks  current  at 
t  started  with  i  customers  in  unit  1  at  t=0  and  r  (<  k  -l) 
service  completions  are  made  before  t.  The  last  service 
completion  before  t  occurs  at  time  u  at  which  there  are  v 
customers  waiting  in  unit  1.  The  number  of  arrivals  between 
times  u  and  t  is  j-v.  Now  there  are  j-1  customers,  excepting 
the  customer  in  service,  at  t  in  unit  1.  Let  the  customer 
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who  is  in  service  at  t  complete  his  service  at  time  v  and 
let  the  service  completion  of  the  j-1  customers  occur  in 
unit  1  at  time  v^.  Finally  we  sum  over  all  choices  of  r, 
v,  u,  v,  vr 

Taking  the  transform  of  (41 )  results  in: 


(42) 


v=0  0  0 


a  H<k)(v) 


♦  V  eW(5)l.i(s)h^k(S)fe-  EV  ^  at 

r=i  o  o 


•  d  H<r)(v) 


r=0  v=l 


o  o 


•  d  H^v)  , 


dt 
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Hence : 

(**3) 


i(b,s,z)  =  E  g^(r)[zh1(s)]vf  e‘svJ 
v=0  '  o  o 

‘  Jt  —fcV. -  h2  (v) 

J  ““V 

r=i  i0  Jo  Jo  j=0 


pb 


•  h*J  (s)  at  d  H*r\v) 


k-1 
+  E 
r=0  v=l 


•  E 
j=v 


— n^rr— 


r  y. 

hg  k  (s)dt  dH^(v) , 


Taking  the  summations  inside  the  integrals  is  justified 
by  LebesgUe  Dominated  Convergence  Theorem.  In  (^3)  to  sum 
the  series  inside  the  integrals  we  use  the  theorem  in  Appendix 
A,  by  taking  y  =  h2(s)  ^  1  for  R(s)  >  0.  Then: 


pflM 
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k-1 

^(S.s.z)  =  |  E  [h^1(s)[umh2(s)>ll[|+X-s-X(Dmzh1(s)h2(s)]] 
m=0  *" 

•  [“X<s>-i5i*4v)(5^vhi(s>h2(s)]V 

v=0 

.  [h2  ( S ) -h2 ( £+X- Xu»mzh1  ( s  )hg( s ) )] 

+  I  g^(l)Ch2Vs)-h2(5+X-Xcinzh1(s)h2(s))] 


+  ^  E  g^)(?)CVh1(s)h2(s)3v[amh2(8)]r 


r=0  v=l 


'[l-h1(e+X-Xa^h1(s)h2(s))|h1(s)]j  ] 


E  rh2"1(s)Ccmh2(s)-l][^X-s-Xumzh1(s)h2(s)]’l 
tt=C  L 

[am[h2 ( s )  - 1  ]  -[g[k  \  §> ( 8  )hg  (  8  )  ) 


•  DigCs )-hg( £+X-Xujzh^(s)h2(s))3 
k-1 

+  E  gf  )(5,0)[h2(B)-h2(^X-Xu)mZh1(8)h2(^))] 
r=i 

+  E  [ch  (s)f  zhi(8)h2(B))-g[r-  ^U,0)  3 

r=0  m  i 

•  [l-h1(e+X-Xu)szh1(s)h2(B)}|  h1(s)3]  ]  » 
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Formula  (26)  gives: 


6^  (5,  wmzh1( 8  )h2(s )  )hg(  S+X-Xiy*^  (s  )h2(  s ) ) 


k-1  r  \ 

+  ^  E  8f;(^0)h2(5+X-Xa)mzh1(3)h2(s)) 


(45)  =  q1(5,umzh1(s)h2(s))  , 

Again,  using  (8b)  we  have: 

g|r^(  5>umzhi(s)h2(s)  )-g|r\  ^,°) 

'  [  “B  2  Vs)  *2<a>I 

f  rh1(r+X-Xcnzh1(s)h2(s))1r'v 


h  (C+X-Xezh  (s)h2(s))-,r 


-  r  r  2 

'  v=i  ri^TjK^rr 


■]  8fJ'('.>°)  , 


Using  this,  the  last  term  in  (44)  becomes  after  simplification: 

(46) 

rJ0  [wmhZ^]  [4r)(5’Vhl(s)h2(s))*gir)(^0)] 

•  [l-h1(s+X-Xomzh1(s)h2(s))  |  h1(s)J 
fh  (s)-h  (E+X-Xc»  zh  (s)h  (s))l 

J  {cvVa)h2(,)} 


!1 

zh. 


k-1 


v=l 


^<mh2^s^V8iV^£'»*0)'«[k^>tJrathi(8)h2^8^h2(8)}  » 


2k 


Substituting  (45)  and  (46)  in  (44)  and  simplifying  we  prove 
the  lemma. 

Lemma  1.7 

For  R(s)  >0,  R(£>)  >  0  and  |z  j  <  1  the  transform 
®i(5>s,z)  is  given  by 

oo  CO 

(47)  ■1(?,«,z)  *  I  r  J“}(0*  (f,i,z) 

n=0  v=l 


=  I  r,  (O'!'  (5jS,z) 

O  IV  V 

v~l 

where  ( %»s>.z)  foi  i  >  1  is  given  in  lemma  1.6,  and 
Qr^  (•)  is  defined  in  (30). 

Proof: 

Upon  taking  transform  in  (39)  we  obtain: 

CO  00  ,  . 

“lj(5's)  *  Jo  Ji  o’iv  (F)  V!'s) 


=  'oriv(S)  *vj(!”3) 
v=l 


Multiplying  both  sides  by  z^  and  summing  with  respect  to  j  we 
get  (47). 

Theorem  1.1 

For  R(s)  >0,  R(i,)  >0  and  \z  |  <  1  the  transform 
^(£>8**)  of  the  joint  distribution  eij(t>x)  of  lueue  length 
and  virtual  waiting  time  at  time  t  for  the  tandem  queue  with 
non-zero  switching  rule  is  gi\en  by: 
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^(Ijs.z)  =  a^^z)  +  i^rxv^)'  [1+x  ®1(^»s>z)]  > 
where  8^(|,s,z)  for  i  >  1  is  given  by  lemma  1.7. 


Proof: 


The  transform  of  (36)  yields: 


(*9)  0**O»s)  =  *~(|,s)  +  »1(5)^,(S.«> 


8oj  I  e~5t  pHk)(t)=°u(k)(oH] «  , 


The  Kronecker  delta  in  the  last  term  is  due  to  the  fact 


that: 


p{l(k)(t)=j,Ti[k)(t)=o|^(oH} 


»  0  if  i  {  0  , 

=  p{^k)(t)=o  |  |(0)*ij  if  J-0  , 

If  M( • )  is  the  renewal  function  of  the  general  renewal  process 
formed  by  the  ends  of  busy  periods  and  m(£,)  its  L.S.T., 
then: 

(50)  j  e_?fc  p{^k)(t)=o  j  ^k)(o)»i}  at 

0 

r00  -ct  r*  -x(t-u)  .  . 

-  J  «  1  e  d  M(u)  » 

o  o 


=  #  • 


where  u  is  the  end  point  of  the  last  busy  period  before  time 
t  and  no  customer  arrives  between  u  and  t. 

Since  the  input  is  Poisson,  between  two  successive  busy 
periods  there  is  a  negative  exponential  idle  period.  Consider 
the  renewal  process  formed  by  the  end  points  of  busy  periods 
and  let  F^*)  be  the  distribution  function  o*  the  initial 
renewal  and  F(*)  be  the  common  distribution  function  of  other 
renewals.  Then: 


F1(x)=G^i^(x),  which  is  the  i-fold  convolution  of  G(0 
F(x)  *  j*  [l-e‘X(x'u)]  d  G(u)  , 


where  G( • )  is  the  distribution  function  of  busy  periods. 
Hence  the  renewal  function  M(t),  which  is  the  expected 
number  of  renewals  in  [0,t],  is  given  by: 

M(t)  =  (F  *  T  F(n))(t)  , 

1  n=0 

Taking  L.S.T.  we  get: 

GO 

n(S)  =  f,U)  I>  f®U) 

n=0 

fl^  ,  R(5)  >  0  , 

=  l^fTCT 


where  fjU)  and  f(?)  are  the  L.S.T.  of  F1(- )  and  F(-) 
respectively,  which  are  given  by  f^(t,)  =  Yi(§)  and 
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5.  Distribution  of  Virtual  Waiting  time 

The  stochastic  behavior  of  the  process 

{T]p^(t),  0  <t  <  “}  is  as  follows:  Tl^k'(0)  is  the  initial 

occupation  time  of  the  server.  If  T|^(0+)  =  0,  then  the 

server  is  idle  at  time  t  =  O  and  until  the  arrival  of  a 

customer  who  initiates  a  busy  period. 

Consider  the  arrival  times  t^,  tg,....  within  a  busy 

period  which  started  at  t  *  0  with  i  >  0  customers  in  unit  1. 

tR  is  the  n-th  arrival  point  of  the  busy  period.  Let 

i  =  mk  +  r(0<r<k,  n  is  a  positive  integer)  and  be 

the  service  time  of  the  n-th  customer  in  unit  v,  v=l,2. 

Then  at  t^  the  arriving  customer  has  a  service  time 

in  unit  1.  Hence  at  t  +0  the  virtual  customer  has 
i+n  n 

to  wait  a  further  X^  units  of  time  more  to  enter  service 

i+n 

in  unit  1,  provided  i+n  is  not  a  multiple  of  k.  That  is,  if 

(k) 

i+n  is  not  a  multiple  of  k  then  at  t^  T|£  '(t)  has  a  jump  of 

Cl) 

Magnitude  X'  .  On  the  other  hand  if  i+n  is  a  multiple  of  k 
i+n 

then  the  virtual  customer  has  to  wait  until  the  completion 

(k) 

of  that  cycle  and  hence  at  t^  \t)  has  a  jump  of 
magnitude: 


d)  +  x(2)  +  x(2 

i+n  Ai+n-k+l  i+n-k+2 


+  X 


(2) 

i+n 


This  is  shown  in  Figure  1. 
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Figure  1 

Graph  of  the  Stochastic  Behavior  of  the  Process  ^[^(t),0<t<5»| 
For  i  >  1,  let  jW^t.x)  =  P-[Tl[k\t)<!c  |^k^(0)=ij  be  the 
distribution  function  of  virtual  waiting  time  P,^(t),  given 
that  at  t=0  there  were  i  customers,  ^(tjs)  is  the  L.S.T. 
of  jW^tjx)  and: 


(53)  ,W**(f;,8)  =  e" ^  ^(t.s)  dt 

o 


Theorem  1. 2 

For  R(s)  >  0  and  R(s)  >  0,  the  transform  (E|,s)  of 
the  distribution  function  jW^(t,x)  of  the  virtual  waiting 
time  T|j^(i)  is  given  by: 


(54) 
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^>1 . 


i  >  1  , 


■where  ^  (f,s)  is  given  by: 
(55) 


1A^#(E,s)  =  |  lJb2"^(s)[ii|ah2(s)-l][f+X-s-\c.^h^(s)h2  (s)]J 
m=  u 

■  [%,D>‘(«)-1]  Kf1(«)f2(3)]1  -  vi(0 

*  %  °r»<e)  3!  (i-«(s)giv)(f’o)}  i 


Proof: 


Prom  Theorem  1.1  the  transform  of  the  joint  distribution 
(^(tjx)  of  queue  length  and  virtual  waiting  time  is  given  by 
©i(c,s,z).  Hence  the  transform  of  the  distribution  of  virtual 
waiting  time  is  obtained  by  taking  the  limit  z  -+1  in 
®^(^>»s,z).  That  is,  from  (48): 


XW.  (f,s)  =  lim  0i(t,s,z) 
z-1 


=  ii(?’*1)  +  FftM1+xa 


1(5,s,l)] 


Hence  it  suffices  io  prove  chat: 


(57)  ^,8,1)  -  1Ai  (r,s) 
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Lemma  1.7  gives: 

(58)  aje,*,!)  =  z  1  oq£j}(0  ij(s,s,i) 


n=0  j=l 


where  fron  lemma  1.6  ^(5,a,l)  is  given  by: 


V^S,1)  =  5  *Q  [h^1^'^mh2(s)-l][§+x-s-Xcinh1(s)h2(8)]] 


[affith2(s)'l}{[«mh1(8)h2(8)]J-qj[|,mnh1(a)h2(s)] 

+  Ji  (l-°m)h2(s)«jV)(^°)}  ] 

Substitution  of  this  in  (58)  leads  to: 
k-1 

V5-8’1*  -  E  jJh2'1(sH»mi>2(s)-i]U*x-5-v.mhl(s)h2(5)3]'1 

{ol[n)C5,^l>1(s)h2(E)]-oq:[n>(6,0) 

•  0^”*l)C?,Ci,hl(s)h2(s)} 

+  A  A  {1"$  h=(8)  bj',)(s’0)}  1 

which  establishes  (57) ,  noting  that  from  (33): 

0^0)(«  =  Bij  8nd  from  ^ > :  2  o<^n'(?,0)  *  Yi(§)  . 
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Limiting  Behavior  of  the  Virtual  Waiting  time  Proce3s 
The  limiting  behavior  of  the  distribution  of  virtual 
waiting  time  is  given  by: 

Theorem  1. 3 

If  X  +  X  a2  <  1  the  limiting  distribution 

lim  jW  (t,x)  =  -W(x)  exists.  The  L.S.T.  of  ^W(x)  is 

t  ® 

given  by: 

tt-ft 

(59)  ^(s)  =  (l-Xc^-Xa2)[l+\  ^  (0,s)3 

where  ^  (0,s)  is  given  by  (55). 

If  Xo^  +  Xa2  >  1,  then  lim  1Wi(t,x)  =  0  for  all  x. 

In  order  to  prove  this  theorem  we  first  show  that: 

Lemma  1.8 

if  p.o(t)=  p[ii[k)(t)  =  0  |  t,(k)(o)  =  i} 
then  the  limit  lim  P.  (t)  =  P*  always  exists.  We  have: 

t  ->co  10  ° 

„  1  -  Xa,  -  Xa„  if  Xa,  +  Xa-  <  1 

(60)  P  =  (  ^  2  “12 

0  if  Xa^  +  Xa2  >  1 

Proof: 

If  m(’)  is  the  renewal  function  of  the  general  renewal 
process  formed  by  the  ends  of  busy  periods,  then: 

picft)=  P^k)(t)  .  0  \  ^><0)  *  i}  -  {  e'X(t_u)  dM(u) 
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By  Smith's  Key  Renewal  Theorem, 


ii®  Pio<L)  =  ^  j  e"AU  du  if  +  ^a2  < 


t  CO 


=  o  if  Xo^  +  Xa2  >  1 

where  >1  =  ^(l-Xc^-Xo:  )'  is  the  raean  renewal  ti®e-  Hence  the 
lemma  is  proved. 

Proof  of  Theorem  1.3; 

Summing  equation  (38)  with  respect  to  j  we  get  the 

distribution  of  virtual  waiting  time  as: 

ft 

(6l)  -jW^tjx)  =  ^.(^x)  +  J  1A1(t-u,x)  dML(u) 


P{^k)(t)=0  |  t>(k)(0)=i}  U(x) 


where : 


(&■)  -.A ^(tjx)  =  E  f  (t,x) 

0*0 

Taking  L.S.T.  of  (6l): 

(63) 

2W*(t,s)  =  xA^(t,s)  +  j  1A^(t-u, s )  dM^u) 


+  p{^k)(t)=o  j  $(k)(0)  =  i} 

Using  Smith's  Key  Renewal  Theorem  (Theorem  U,  Appendix  D) 
and  lemma  1.8  and  taking  the  limit  of  (63)  we  have: 


(64) 
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lim  ^W*(t ,  s )  *  X(l-Xa1-Xa2)|  ^(ujsJdu+^-Xa^-Xag)  , 

t  -V  ®  o 

if  1-Xa^  -  Xa2  >  0 
*  0  if  1-Xa^  -  Xa2  <  0. 

00 

*  -  V 

since  from  (4l)  it  can  be  shown  that  lim  T  .(t,s)  *  0 

t  »j=0 

*  00  * 

which  implies  from  (39)  that  lim  -A.(t,s)  =  lim  E5.(t,s)=0. 

t  o  1  t  j=0 

Again,  (64)  can  be  written  as: 

(65) 

-X* 

lim  ^(tjs)  =  (l-Xa1-Xa2)[l+X-1A1(0,s)]  if  l-xa^-Aag  >  0, 

t  -'CO 

=  0  if  l-XOj-Xaig  <  o 

From  (55)  it  can  be  shown  that  Ax  (0,s)  is  continuous  at  s=0. 

Hence  by  Zygmund's  Theorem  (Theorem  1,  Appendix  D)  the 
limiting  distribution  lim  ^W(t,x)  =  ^W(x)  exists  and  the 

t  cn 

L.S.T.  of  xW(x)  is  given  by  (65). 

Formula  (55)  can  be  rewritten  as: 
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*  I  „r  (0)  £  (l-t£)h*(s)g(v>(0,0)}  1  ) 

j=l  v=j  ra  c  0  J  J  7 


=  lX(s)  +  ’J2(s) 

where  ^(s)  and  ;’2(s)  "re  respectively  the  first  and  the 
second  term  of  (66) 

Taking  the  limit  s  0  in  (66),  we  see  that  the 
numerator  of  £2(s)  is  zero  while  its  denominator  is  non¬ 
zero.  Hence: 

jAj*(0,0)  =  ^(0+) 

=  °j  +  a2 

1-Xa^  -  Xa2  * 

which  together  with  (59)  gives  ^c:(0+)  =  1. 


Expected  Value  of  the  Limiting 
Distribution  of  Virtual  Waiting  time 


Let  denote  the  expected  value  of  the  limiting 

distribution  of  'O^(t).  Then 

(67)  M^Ck)  ■  -  fj 

■  -  »<i-*vxa2)  h  iAr(°’s>]s=0 

From  (66): 


where  the  number  of  primes  indicates  the  number  of  successive 
derivatives  taken  with  respect  to  s.  Let  A^'s)  denote  the 
numerator  and  A2(s)  the  denominator  of  ^(s),  so  that: 


Applying  de  l'Hopital's  rule  four  times  on  the  right  hand 
side  of  (69)  we  get: 

jVc*)  -  [i”(0)4"'(0)  -  ^(0)4-  (0)]  [i"(0)J 


where  after  simplifying  we  obtain: 
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A*(°)  =  +  o.2) 

^(°)  =  "3k[«2^l+ ^^ia2^2^+^k"1^a2^1+0f2^+^2^al+0!2^] 

^(0)  =  "  2koi2(-l  |-  X  +  Xoig) 

A2(°)  =  3kj^2(k-l)o;2(-l+Xa1+Xo:2)+3Xa2(01+3D!1a2+32) 


+  ^gC-l+Xo^+XagjJ 


Hence: 


(70) 

1  2(7oraj7tojr 

Differentiating  ^(s)  with  respect  to  s  and  setting  s=Os 
a?  k-l  cv  . 

(0)  =  -f  T  - U-l 

2  X  m=l  (l-a  )2  1  m 

v  nr 


+  ?,  0ri j(°)  V  f1’#  b5v)(o,o)}  , 

J=1  V=0 


(71) 


-a0  k-l  a' 

_ £  r  m 

X  «-i  1-® 

m=i  m 


K-!  K-l  >)(0>p) 


a?  k-l 

*  T  1;1  orij(0)  -  /  (1:7^ 

J-l  v=j  m»l  v  m 


From  the  properties  of  tht  roots  of  the  equation  z  -1  =  0 


we  find  that: 
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k-1 

L 

m=l 


l-ia_ 


m 


k-l 

2 


(72) 


k-l  a  k-l 
•;  m_ 

‘  1-cu  -  . 

m=l  m  m=l 


m 


By  the  method  of  partial  fractions  we  get: 

a£+1  s  l-(v+l)(l-a)m)+v(l-<um)2+(v-l)a)ra(l-a)m)2.(v-2)c^(l-om)2 


_  v-2/.  v2  v- 1/ ,  \2 

h  2uJm  (l"u)m)  +  am  (l'uk) 


which  gives: 

k-l  s)  (1-ul) 

m  m 

m=l  (l-u)m)2 


k-l 


k-l 


k-l  a) 


v+1 


m 


m=l  ^”cm  m=l  m=l  (l-u  Y 


m  k-x  r 

m=l  L1" 


- v-(v-l)ui  -(v-2)a' 

id  v  '  m  '  '  m 

m 


-  2a 


v-2  v-11 

j  -a 

m  m  J 


v(~)-v(k-l)+[(v-l)+(v-2)+  •••■  +  2+1] 


k-l 

r 

since  a  - 
m=l 


1  for  1  <  r  <  k-l. 


That  is: 

,  .  O  (l-CJ^)  /,  \ 

k-l  nr  m  v(k-y) 

(73)  ^  ' '  2 


Substitution  of  (72)  and  (73)  in  (71)  yields: 
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(7M 


,  (K-D  ct2 
>jj(°)  -  —nr 


a2k-i 


k-1 


(v), 


•jr  I  r  (0)  2  v(k"v)  s:  (0,0) 

Zx  j=l  0  13  v=j  J 


Formula  (68)  together  with  (70)  and  (74)  leads  to: 
(75) 


a  **,  x  1 

^  iAi  (0>s)  = 


2(l-\a^-Xa2) 


£(k-l)o:2(al+a2)(l“\0!1-Xa2) 


(K-l)o, 


2X 


a  k-1  k-1  ,  . 

of  '  or!,(0)  v(k-v)g  v)(0,0) 

a  j=l  0  lj  v=j  J 


It  follows  from  (67)  that: 

X(01+^1a2+02)  (K~l)arz 

k-1  = 

'1 


(76)  \W  -  2(i;^2)- 


a, 


k-1 


k-1 


(v>, 


+  -f  (l-Xc^-Xag)  ^  ^(0)  T  v(k-v)gj  '(0,0) 


where  Qr^.(*)  Is  defined  in  (30). 


It  is  worth  noting  that  for  k-1: 

MT^(1)  =  ^^Bi+2aia2+^2^  I  2(l-Xa1-Xa2) 


which  is  the  expected  value  of  the  limiting  distribution  of  the 
virtual  waiting  time  of  an  m|g|i  queue  with  impute  rate  X  and 
service  time  distribution  H,*H2(‘)-  For  k-1  the  tandem  model 
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reduces  to  an  M|G  (1  queue  with  service  time  distribution 

VH2(*)  ' 

Conjecture: 

Por  all  k  >  1: 

(77>  < *\(k) 


The  proof  for  k=l  is  simple: 
Prom  (76)  we  have: 


\(2)  *  \(1)  "  I  [^V^ll*^! (0'0)-l] 

<|[(1-V>l-Xa2)orii(°)-l]<0 

03 

3  of  Appendix  B  -^(0)  < 


since  from  Theorem 


6  .  Queue  length  Process 

Let  P. .(t)  denote  the  probability  that  at  t,  j  customers 
are  in  unit  1,  given  that  the  service  started  at  time  t=0 
with  i  customers.  That  is: 

(/8)  P..(t)  =  P^(k)(t)  =  j  j  S(k)(o)  =  i} 

Let  ix.  .(£.)  be  the  Laplace  transform  of  P.  .(t)  and: 
ijv  *  lj 

00 

(79)  «4U,2)  =  ,  |z]<l,R(c)>0 

j -O  1J 

0.0  |z  |  <  1,  R(?)  >  0 

Por  i=o  and  ,j  >  0  we  have: 


(80) 


j  and 
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noi^  =  XH  +  X+f-  *13^ 


(81) 


=  ikl1  +  X  *lU>l)] 


The  equation  (80)  is  obtained  by  considering: 

(i)  If  j=0  then  there  can  be  either  no  arrival  in  [0,t] 
or  there  is  a  negative  exponential  idle  period  followed  by  a 
busy  period. 

(ii)  If  j  >  0  there  is  a  negative  exponential  idle  period 
followed  by  a  busy  period,  and  in  this  case  the  first  term  on 
the  right  hand  side  of  (80)  vanishes. 

For  i  >  0  we  have: 

Theorem  1.4: 

The  generating  function  ( t, , z )  is  given  by: 

(82)  *  xt(5.»)  *  ■$£-*(%)  [L  *  X  xl(C’z)] 

where 

r  -I'l 

(83)  Xi(?»z)  =  |^(  >j+X-Xz '  C?-h^(  E,  -  \-Xz)]  J 

J^s1[z-h1(P>+X-Xz)]-zLl-h1(t,+X-Xz)]  yi(  ?) 

a  k"  1 

-  (z-l)h1(^X-Xz)^or.(§,z)'  I  ^jtUCg^ (?»*)♦  r  8.V^.0)l}  ]» 

j “»  v~l 

and  r.  . ( L,) ,  r.(f,z)  are  defines!  in  (30)  and  (31)« 

o  ij  0  1  - 


kz 


Proof: 

Tha  generating  function  it  (5»z)  is  obtained  by  taking  the 
limit  as  s  -t0+  in  the  transform  6L(e;,s,z)  of  the  joint 
distribution  0^(t,x)  of  queue  length  and  virtual  waiting  time. 
From  Theorem  1.1  xe  have: 

n  (c,z)  =  lim  e  (rD, a, z) 
s  Of 

(84)  =  [w  ^(5,0, z)] 


Kence  it  suffices  to  prove  that: 


(85)  2i(^J0,z)  =  X.U,s) 

From  lemma  1.7,  8.  (-,*,•)  is  given  by: 

(86)  3i.U,0,2)  =  I  V  q(^'r)  V5’°,z) 

n=0  3=1  3 

where  0,z)  from  lemma  1.6  is  given  by: 

J 


(87) 


,U,0,z)  =  |  (5+X-Xz)[z-h1(§+X-Xz)]] 


[Cl-h1(  %+\- *z )  lz  j+1-  [z-h1(  |:-a-Xz  )]q^(^,z) 

k-1 

+  '■!-l)h1(5+X-Xz)^^(§,z)+  7  g^(§,0)j^ 

since  other  terms  iu  the  -nidation  on  the  right  hand  side  of 

(40)  for  m=l,2, - ,k-l  vanish  as  s  -■  o.  Substituting  (87) 

in  (86)  we  find  that: 


'+3 


(88)  ^(5,0, z)  »  jjs+X-XzKz-h^EtX-Xz)]] 

r.Cl-h^E+X-X.)]  £  lo4n;(5,2)-olin)(5,0)] 

L  n=0 

-  [z-h  (s+x-x>)]  r  0<4”+1>((,,z) 

n=0 

+  (z-l)h  (t;+x-Xz)  r  r  (dg,W(5,zK  I  g^v)(l,0)}  1 
j=l  v=j  J  J 

=  [(^+X-Xz)[z-h1U+X-Xz)]l 

|"z1[z-h1(5+X-Xz)]-z[l-h1(C+X-Xz)]  Yi(0 


(z-l)h1(|+X-Xz)  |or.(^,z)-  I  ori;JU)  [gjk^(5,z) 

g(,v)(5,o)]}  ] 


k-1 
+  £ 
v~l 


which  proves  (85). 

Limiting  Behavior  of  Queue  length 

The  limit  of  P.  .(t)  as  t  -•  <=>  always  exists  by  a  thearem 
1 J 

of  Smith  (1955). 


Let  us  denote: 


(89) 

and 


p .  =  lim  P.  ,(t) 
J  t  o  lj 


(90) 


P%) 


j=0 


Pj 


1*1  <  1 


> 
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Then  by  a  standard  Tauberian  Theorem  (Theorem  5>  Appendix  D) 
we  have: 

(91)  P*  =  lin  M  P..(t)dt 

3  t  t  3o  1J 

CO 

=  lim  5  [  e-^  P  .(t)  dt 
5  -j  0  o  iJ 

which  together  with  Theorem  1.4  gives: 

P*(z)  =  lim  5  n ,(F.,z) 

§  -t  0 

*  ?Um0^tT  tw 

(l-Xa,“Xa2)[l+XX1(0,z)]  if  l-Xa^-Xag  >  0 

(92)  ={ 

0  if  l-Xa^-Xctg  <  0 

where  x^(0 ,z)  from  (83)  is  given  by: 

<93>  V°'’>  -  {-^1(V-X.{orl(0,z) 

-  I  /..(oiufio,!).  T  f5v)(0.c)4]  } 

j=l  J  J  V=1  J 

Substitution  of  (93)  in  (92)  yields  the  following: 

Theorem  1.5 

*.  . 

The  generating  function  P  (z)  of  the  limiting  probabilities 
*  .  . 

P.  of  P. .(t)  is  given  by: 
j  ij 
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P*(z)  =  (l-Xa1-Xa2)h1(x-Xz)[z-h1(x-Xz)]”1  j-1  +  or1(0,z) 

- [«!k)(°,z)+S  e^,(o’o)] } 

J 

where  Qr^(*)  is  defined  by  (30). 

The  Steady  State  Expected  Queue  length 

Let  M^(k)  denote  the  mean  of  the  limiting  distribution 

(s) 

of  I  (t)  as  t  ®.  Theorem  1.5  gives: 
m500  -  fj  p’U)]^ 


(l-^-Xctg) 

2(l-Xa~)S 


CO 

[C2Xce1(l’-Xct1)-8-X2p1D  I  S  {^(0) 


Viv<°> 


*  (1^)n(J-i){0ry(°)-ji0rlv(o)^'(°)}  ] 

(1-Xcc,  “Xo:?)  00  r  f  2  4 

- - — y~  r  |o,lr.2xa1(l-Xa1)^X<ip1Hj-l)(l-XaL)j 

j=iL 


I  r.  .(0)  -  T  r.  (0)  g^k 
lo  lj'  ^  lv'  vj 


(0)  6^(0)}  ] 
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CHAPTER  XI 

SINGLE  SERVER  TAJffiEM  QUEUE 
WITH  ZERO  SWITCHING 

Here  we  consider  a  tandem  queue  with  a  zero  switching 

t 

rule  in  units  1  and  2.  At  t=0  the  server  starts  in  unit  1 

y 

and  continues  to  serve  there  until  the  queue  in  unit  1  be¬ 
comes  empty.  After  completing  the  1-task  the  server  switches 
to  unit  2,  serves  all  the  customers  there,  then  switches  back 
to  unit  1  and  continues  in  this  manner.  If  the  whole  system 
is  empty  the  server  waits  in  unit  1  for  the  arrival  of  a 
customer  who  initiates  a  busy  period. 

The  analysis  of  the  queue  with  zero  switching  case  is 
easier  than  non-zero  switching  case.  The  distribution  of 
busy  period  for  the  zero  switching  case  is  the  same  as  in 
non-zero  switching  case  (Chap.  I  section  2). 

1.  Transition  Probabilities  of  the  Basic 
Imbedded  Semi -Markov  Process 


We  use  the  same  notations  as  in  Chap.  I. 


^7 


(1)  Vx)  ■  Kv**  i  Sh-i  - 4} 


■  rf  f  <)(u)e-X(v'U>  dH<r)(v-u) 

r®i  o  u  * 

(u)(v) 

The  transforms  ^(s)  and  ^(sjz)  of  ( • )  are  given  by: 

(r) ,  .  r®  -tB±\  w 


(2) 


W>  *  £  C<*>  I  *-(stx)x  ^  »'r)(x> 


<l(s,z)  =  r  g(r^(s,0)  hJVs+x-Xz) 
r=i  2 

(3)  *  Y*js,  h2(s  +  X  -  Xz)} 

(by  lemma  1.2) 

(n) 

If  oQij  '*)  are  the  tabo°  Probabilities  defined  in  (1.28)  then: 
-  6ij  u<x> 


and 


00 


Qi.l'(x)  =  r-  fi‘Qiv"1^x"u)d^(u)»  n>l 


4(n) 


v=l  o 


vr 


Their  transforms  are  given  by: 


U 


Af^ *  jx  A?l)u  av^s)  * n  > 1 » 


(5) 
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(a,z) 


i 

z 


(6) 


*  I  q^Jj  1^(a)  q  (s,z),  n  3>  1, 

01  v*l  v 

A  further  discussion  of  these  transition  probabilities  is 

given  in  Appendix  B. 


2.  The  Joint  Distribution  of  Queue  length 
and  Virtual  Waiting  time 

For  the  definitions  of  queue  length  and  virtual  waiting 
tine  we  refer  to  Chap.  I  sec.  4.  Let  §(t)  and  ^(t)  respectively 
denote  the  queue  length  end  virtual  waiting  time  at  t.  Define: 


(7)  ei;J(t,x)  -  p[i(t )=j,  ^(t)  <  x  |  5(0)-i] 

(8) 

SU(t,x)  “  0  <  \(t)  <x,\(i)fO  for  all  T<E(0,t] 

I  1(0)  *  i  } 


Analogous  to  equation  (1.38)  we  obtain: 

ft 

(9)  S1;J(t,x)  *  *i;J(t,x)  +  J  *y(t-u,x)  dM^u) 


+  i{5(t)-J,  \(tbO  |  |(0).i}  U(x)  , 


Let  *  (t,x)  denote  the  probability  that  at  time  t  the  original 

cycle  has  not  ended  that  3(t)  =  j ,  0  <  1\^(t)  <  x  and 
T^(t)  ^  0  for  all  t€  (0,t],  given  that  at  t=0  the  service 
started  in  unit  1  with  i  customers. 

H  j|  || 

We  define  the  transforms  ©^(^s),  $i;j(S,s),  ^(§,8), 
0i(C,s, z),  ^(1,8,2!)  and  as  in  Chap.  I 

Lemma  2.1 

Por  R(s)  >  0,  r(^)  >  0  and  |z|  <1,  the  transform 
given  by: 

(10)  *!(«,»,*)  *  ^x.s.uhi(s)  {«1(5.iT:£)-q1(l.^1(»)) 


z[(lh1(«))1-v‘;5)Jh1(s)-h1(E.X-Wi1(8))] 

^zh1(s)-h1(|+X-Xzb1(s))]  j  ,i>l, 


Proof: 


The  probability  ^  .(t,x)  is  given  by: 
^  J 


(a) 

♦l3(t,x)  . 


if  rr^'w^1  ^ 

r=i  o  t  v 

(u)(v)  (v.)  ,  v  m 

dHg  (v-u)dH^  (v^-v) 

”  ( CC  ^ 


(u)(v)  (vx) 
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The  first  term  is  obtained  by  assuming  that  the  server  is 
performing  a  2-task  at  t.  The  cycle  of  tasks  in  which  the 
server  is  serving  at  t  stores  with  i  customers  in  unit  1  at 
t=0  and  the  unit  1  becomes  empty  after  r  pervices  at  time 
u(0  <u  <t).  There  are  j  arrivals  in  (u,t).  The  service  of 
these  j  customers  starts  after  the  completion  of  the  2-task 
in  progress.  Let  this  2-task  end  at  time  v(t  <  v  <  t  +  x). 

The  service  completion  in  unit  1  of  the  j  customers  occur 
at  tiro  v^(v  <  v^  <  t  +  x).  Now  we  integrate  and  sum  over 
all  choices  of  r,  u,  v,  ar.i  v. . 

The  cocoad  term  is  obtained  by  assuming  that  the  server 
is  performing  a  I-tc.sk  at  t.  The  cycle  of  tasks  in  which  the 
server  is  serving  at  t  starts  with  i  customers  in  unit  1 
at  t=o.  Let  there  be  r  service  completions  in  unit  1  before 
t.  The  last  of  these  occurs  at  time  u  and  at  this  time  there 
are  v  customers  waiting  in  unit  1.  There  are  j-v  arrivals 
in  unit  1  in  the  interval  (u,t)  so  that  at  t  there  are  j-1 
customers  in  unit  1  excepting  the  customer  in  service.  The 
service  completion,  in  unit  1,  of  the  customer  in  service 
occurs  at  time  v.  The  service  completion,  in  unit  1,  of  the 
j -1  cu-J  om.ers  occurs  at  time  v^.  Finally  we  sura  over  all 
choices  of  r,  v,  u,  v,  and  v  . 

Equation  (11 )  is  valid  for  ail  j  >0.  For  j=0  the  last 
term  clirsppcvs.  Upon  t-'k'rg  transform  in  (11)  we  find: 
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(12) 

** 

'id 


g^)(S)hJ(s)jVsvjV'^-=H  ,t  aH(r)(v) 


r=0  v-1 


o  o 


dt  dH^v) 


Hence : 


^(SjSjz)  =  £  L  .(£,s)  zJ 

j-0  1J 


*  Fx-'s-th^s) 


n  co 


+  L0  ^g^)(E,)h^1(s)Z'J[h1(8)-h1(|+X-Xzh1(8))]} 

Using  (3)  and  (1.1*0  and  simplifying  we  prove  the  above  lemma. 
Lemma  2.2 

For  R(s )  >0,  R(§)  >0  and  \z  I  <1  the  transform 
5i(5>s,z)  is  given  by: 

(13) 

sAs-*>  -  vr-s-^c.')  ^x[°4n)  ( s-*T*>v4n) <  ^V3»] 


z[^1  ( s )  -h1  ( E>X-  Xzh1  ( s ) ) 

[zh^ ( s )  -  h1(|+X-Xzh1(s))^j"1 


(n), 
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Proof: 

As  in  lemma  1.7  we  have: 

(14)  $  (*,, a,z)  =  I  I  qj^U)  ♦(&,«,*) 

n=0  v=l 

Substituting  ^(£,s,z)  from  lemma  2.1  and  simplifying  the 
result  follows.  The  convergence  of  the  serices 

I  '(s,z)  is  discussed  in  Theorem  3  of  Appendix  B. 

n=l  01 

Theorem  2.1 

For  R(s)  >  0  R(s)  >  0  and  Jz  J  <  1,  the  transform 

of  the  joint  distribution  8^(t,x)  of  the  queue 
length  and  virtual  waiting  time  at  t  for  the  tandem  queue 
with  zero  switching  is  given  by: 

(15) 

*,*)  =  [l+X^d,*,*)]  ,  i  >  l  , 

where  ^(SjSjz)  is  given  by  lemma  2.2. 

Proof: 

Similar  to  the  proof  of  Theorem  1.1. 

3.  Distribution  of  Virtual  Waiting  time 

The  stochastic  behavior  of  the  process  (T^(t),  0<t<“  } 

my  be  described  as  follows:  T^(0)  is  the  initial 

occupation  time  of  the  server.  If  T^(oO  -  0  then  the  server 

is  idle  at  t=0*-.  Let  i  be  the  initial  queue  length  at  t=0 

and  t  the  n-th  arrival  point  and  the  service  time,  in 
n  l+n 
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unit  v,  of  the  '.astomer  arriving  at  time  t  .  At  t  the  value  of 

n  n 


^(t)  has  a  jump  of  magnitude  X^.  Between  any  two  arrivals 
T^(t)  decreases  linearly  with  slope  -1.  As  soon  as  T^(fc)  reaches 
zero,  it  jumps  suddenly  to  a  magnitude  equal  to  the  tcfcal  str/ice 
time  of  all  the  customers  present  in  unit  2  at  that  time, 
after  which  it  linearly  decreases  with  slop  -1  until  the  arrival 
of  a  new  customer.  This  is  shown  graphically  in  Figure  2. 


Figure  2 


Graph  of  the  Stochastic  Behavior  of  the 
Process  0  <  t  <  ®} 
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Let  ^(t.x)  =  P(nL(t)  <x  |  |>(0)  =  1}  be  the 
distribution  function  of  the  virtual  waiting  time  T^(t), 
given  that  at  t=0  there  are  i  >  1  customers  in  unit  1.  Let 
,.jf*(§,s)  be  the  transform  of  ^(t.x)  defined  in  (1.53)* 
Theorem  2.2 

j| 

For  R(s)  >  0  and  R(£)  >  0  the  transform  (£,s) 
of  the  discribution  function  ^W^(t,x)  of  the  virtual  waiting 
time  (t )  is  given  by: 

(16)  jW^g.s)  =  [l*l 

where 

(!7)  -  [^-s-Xh1(s)]  1{h^(s)-Yi(£i) 


+ 


5+1-s 

X 


)-4(s. 


X 


and 

aQ(l,z)  =  z 

(18)  anU,z)  =  ^{i,hgt|+X-Xa^_^( 5>z)]}  ,  n  >  1, 

Proof: 

We  get  ^^*(£,8)  by  taking  the  limit  z  -c  1  in  the  trans¬ 
form  0^5, s,z)  of  the  joint  distribution  t,x)  of  the 
queue  length  and  virtual  waiting  time.  Theorem  2.1  leads  to: 


s)  -  lim  91(i,s,s) 


z  ■*  1 


(19) 


si(fe,a,1)  *  5^^fef[1+X  *1(?,8,1)] 


Hence  it  suffices  to  prove  that: 


w  y 

(20)  fii(§,s,l)  -  (S,a)  ,  i>l. 

From  lemma  2.2  we  have: 


(21) 

^(6.8,1)  -  U+X-a-Xh1(s)3‘1  |h*(s) 
-  t  04n)(5.\(?))} 

n=o  ■* 

Equation  (3)  gives: 

o^tS-VO)  .o9^n)(S,0)  + 

•  0«iB)(t.o)  * 

■  0«ia)(S.O)  + 


n=l  0  1  K 


a5, 


4ntl\l,  ¥■) 


which  upon  summing  over  n  leads  to: 

r  <4")<s.Vt)) '  *  0<4n)(s.°)  *  r0q<n)(?,^) 

n=o  •L  n=o  °  a  n»l°  *  K 

■  1,1  (?)  *  T  01<n,t5.^) 

n»l 

(For  ccnvergence  of  these  jeries  we  refer  to  Theorem  3  in 
Appendix  B) 

Substitution  of  this  in  (21)  yields: 
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Using  "lemma  2  of  Appendix  B  and  simplifying  the  above  expression 
we  prove  (20), 


Limiting  Distribution  of  Virtual  Waiting  time 
Let  ^W(x)  be  the  limiting  value  of  ^(t,x)  as  t  -t  ® 
and  let  ^ia(s)  be  its  L.S.T.  The  existence  of  the  limiting 
distribution  can  be  proved  by  Zygmund's  theorem  as  in 
Chap.  I  section  5« 


Theorem  2.3 

The  L.S.T.  of  the  limiting  distribution  ^W(x)  of  the 
virtual  waiting  time  T^t)  is  given  by: 

(23) 

!<■>  *  i  [w>.¥)]  i 


if  l-\o^-Xa2  >  0 


=»  0  otherwise 

where  the  functions  an(,,>)  are  defined  in  (19). 

Proof: 

Similar  to  the  proof  of  Theorem  1.3.  Aa  in  (1.65)  we 


obtain: 
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(24) 

]_o( s )  *  (l-Xo^-XagJtl  +X  1A1  (0,s)]  if  l-Xc^-tag  >  0, 

=  0  otherwise 

where  ^  (0,s)  is  given  by  (17): 

,C(o.«)  *  cx-s-M^is)]-1  k:s)-i  +  t  u„(o,iji) 

-1 

-an(0,l)]} 

<25>  *  -s-^kTJT  W‘> +  * 

i  n=x 

The  convergence  of  the  series  £  [l-a  (0,~£)]  is  discussed 

n*l  n  K 

in  lemma  7  of  Appendix  B. 

Substituion  of  (25)  in  (24)  proves  the  theorem. 

Taking  the  limit  as  s  ■*  0 in  (23)  we  observe  that: 


=  1  if  l-Xa^-Xctg  >  0, 


by  lemma  1  of  Appendix  B. 

Expected  Value  of  the  Limiting  Distribution 
of  Virtual  Waiting  time 

Let  denote  the  expected  value  of  the  limiting 
distribution  of  "^(t).  Taking  the  derivative  of  (23)  results 
in: 
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\ 


d  Lu(s) 


ds 


1 


8-0 


(l-Xo'-XaJ  f  ® 

(27) - 7k  *  lx  p,  [1  +  I  a'(0,l)3 

r  X  x  n=l  n 


2(1-^) 


^  \  <(«>.«} 


Using  lemma  1  of  Appendix  B  and  simplifying  (27)  we  obtain: 

x  fc1+02+ 2  ^r^r^aia2  ^ 

(28)  M-  =  -  ^ 


TL  \Qn  2 


i 


Computation  of  higher  moments  seems  to  bu  very  tedious. 
X(31+02+2oLa2) 

If  we  denote  =  2(TT — I XcT) —  wtl^ch  *s  ^he  steady  state 
expected  vjrtual  waiting  time  of  an  M  |G  |l  queue  with  service 
time  distribution  H1^iHg(*)»  then  it  can  bo  shown  that: 

1  °i  11 

2^  -  1  <h\  <Cl  \  "  C2  °i>  for  “Vc1>2»c2<2» 

U.  Distribution  of  Queue  length 
Let  Pi;j(t)=P{§(t)»J  |  |(0)=i}  and  n^?)  be  its  Laplace 


transform  and 
(29) 


^(5,*)  -  r  ,  b(0  >  o,  I*  |  <i 


j=0 


or  3(5)  >0,  |tl  <  1 


Theorem  2.1* 


The  generating  function  «^(?,i)  is  given  by: 


■  xi(6>2)  *  Cl  +  i  VH»  ,  t  > 


where  X^(£,z)  g£ven  by. 

(31) 


♦  z[*-h1(^x-U)]'Yl-hJ,(5.x-!ls)] 

Jo  [*»(?,I)  '  a»(5'V«»]  }  .1>1. 

and  the  functions  aj-,.)  are  defined  in  (18). 


Proof: 


Analogous  to  equation  (1.84)  we  obtain: 

(32)  *  {i(5’0ll)  +  5 *i(5,0.» 

where  Sj^O.z)  is  obtained  from  lemma  2.2: 

+  *[*-h(5+»-»«)]'1[l-hI(5*x-u)] 

4W*>(t->  •  A" C!.\C*)>]  } 

*  \(?,x) 

(by  lemma  2  in  Appendix  b) 
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which  together  with  (32)  proves  (30). 

Limiting  Distribution  of  Queue  length 
Let  P*  =  lim  P,,(t) 

J  t  -{•  OD  1>J 

and 

(33)  p*(z)  =  r  P*  ,  |zj  <1, 

j*o  0 

where  the  existence  of  the  limit  is  established  as  in  Chap.  I 
section  6. 

Theorem  2.$ 

The  generating  function  P  (z)  is  given  by: 

(34) 

*  (1-Xa.  -Xa9)h  (X-Xz)  r 

p  <»>  ■  -iffhsFr —  J0 11-^0.-)]  ^  > 0  > 

»  0  otherwise , 

where  the  IVacfcicns  • , • )  are  defined  in  (18). 

Proof: 

It  follows  as  in  (1.92)  that: 

(35)  P*(*)*  lim  r 

»  (l-Xc^-Xag)  [l  +  X  VO.*)]  if  1-X^-XOg  >  •, 

■  0  otherwise, 

where  V°»*)  is  given  by  (31): 


61 


(36,  X^o,.)  =  [.nC0,i)  -  an(o,z)] 

+  2[z-h1(x-\z)]  [l-h^X-Xz)] 
00 

J0  h(0’z)  -  V0-1’] } 


-1  r  MX-Xz)  » 

■  7  ll  +  ^5^537  JQ  [W°*z>]  } 

Substitution  of  (36)  in  (35)  proves  the  theorem. 

It  can  be  shown  that; 

P  (1)  =  1 

Ebcpected  Value  of  the  Limiting  Distribution 
of  Queue  length 

Let  denote  the  expected  value  of  the  limiting  distribution 
§(t).  From  Theorem  2.5  we  obtain: 

M? '  h 

(l-X^-Xag)  ,  « 

■  Si]  Jo  <i0A) 

*  a-^)  }  •*(0,1)} 


X  (&^P2)^XOi(l-Xai) 


(by  lemma  1  in  Appendix  B) 
■  X  Oi  +  X 


(37) 
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5.  Distribution  of  Total  Time  in  the  System 
Consider  all  the  customers  present  in  unit  1  at  time  t. 
Their  total  time  required  to  complete  services  in  unit  1  as 
veil  as  in  unit  2  is  defined  as  the  total  time  in  the  system 
at  t.  We  denote  by  ^(t)  the  total  time  in  the  system  at  t 
for  the  sero  switching  case. 

Let  Wt( t,x)  =  P{T|2(t)  <  x  |  §(0)  =  i}  be  the 
distribution  function  of  the  total  time  in  the  system  T)2(t), 
given  that  the  queue  length  at  t=0  is  i. 

Further  ve  denote: 

(38) 

gA^tjjO-F^O  <  ^(t)  <  x,T12(t)/0  for  all  T6(0,t]  \  |(0)  =  ij 

and  let  2X^(t,x)  be  the  probability  that  at  time  t  the  original 
cycle  has  not  ended  and  that  0  <  T]g(t)  <  x  and 
^(t)  i  0  for  all  T€(0,t],  given  that  at  t*0  the  service  started 
in  unit  1  with  i  customers.  Analogous  to  equations  (1.38)  and 
(1.39)  ve  obtain: 

(39) 

ft 

2Wi(t,x)  -  ^(t.x)  +  J  gA1(t-u,x)  dM1(u) 

0 

♦  p  (\(t)  -  o  \  5(0)  -  i]  u(x) 

and 
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For  R(s)  >0  and  R(|)  >0  we  define  the  following  transforms: 

#(W>  -  £  I/“  «*<*•*> 4t 

2a“(M)  ■  i  «'9’  Joe-“  « 

*-«..>  -  r  c  e"x  d*<t>x)  at 


Lemma  2.3 

For  R(s)  >  0  and  R(§)  >  0  the  transform  ,  (§,s)  is 

"  o 

given  by: 

(41) 

nlf  r-g+XY  (s) 

jZ*U,b)  =  [e-s+xY^sJ-Xh^sJY^s)  J*|qjU» - - ”3 

-  qjC^.hgCsjY^s)]  + 

--1 

[Yi(s)-hi(!+X-Xh2(s)Yi(i))J 

[h1(s+X-XY1(s))-h1(5+X-Xh2(s)Y1(s))] 

[(h^v.))* .  4s.zp)] } .  j  >  i  • 

Proof: 

The  probability  ?X^(t,x)  is  given  in  ttrms  of  the 
probabilities  G^(0  by: 
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(42) 


,t+x„t+x  „t+x 


r  r  r  r  f  f  [  [  a  g  */(u) 

3  r~j  v, aOv9=Or. =v. +v_  o  Jt  v  J0 

1211  2(u)(v)  (vL)  (Vg) 

.  e“X(t-u)  ftfcu)!  j  e-X(v-t)  Lx(y:t)l  2  d  H(r)(v.u) 


V 


2  * 


(r,)  M 

•  dv1QVv2,o(VT>  dv  H2  <Vt1> 


00  CD 

+  I  F 


0C 

I 


00 

F 


GO 

F 


-t  „t+x  -t+x  «t+x 


r*0v^lv^=0v230r^=v4-v^+Vg-l  o 


rr 


U 

O  u  V  V. 

(u)(v)  (VjHvJ) 


*s)<,> 


(ri)  (r+v+v,)/  \ 

•  d  G  J-  .  ,  (v.-v)  d  H„  1  1  2  V 

Y1  v+v]_^v2_^»°V  1  v2  2 


The  first  term  is  obtained  by  assuming  that  the  server  is 
performing  a  2- task  at  time  t.  The  cycle  of  tasks  in  which 
the  server  is  serving  at  t  starts  with  j  customers  in  unit  1 
at  t«0,  and  the  unit  1.  becomes  empty  after  r  services  at  time 
u.  The  number  of  arrivals  in  unit  1  between  times  u  and  t  is 

The  r  customers  in  unit  2  at  time  u  have  service  completion 
at  time  v  and  v2  is  the  number  of  arrivals  in  unit  1  between 
tires  t  and  v.  At  time  v  there  are  customers  in  unit  1. 
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Starting  with  v^+v2  customers  at  time  v,  unit  1  becomes  empty 
after  r^  services  at  time  v^.  The  service  completion  of  the 
customers  in  unit  2  occurs  at  time  Vg.  Finally  we  integrate 
and  sum  over  all  choices  of  r,  v^,  Vg,  r^,  u,  v,  v^,  and  v^. 

The  second  term  is  attained  by  assuming  that  the  server 
is  performing  a  1-task  at  time  t.  The  cycle  of  tasks  in  which 
the  server  is  serving  at  t  starts  with  j  customers  in  unit  1 
at  t»0.  There  are  r  service  completions  in  unit  1  before  t. 

The  last  service  completion  before  t  occurs  at  time  u,  at 
which  there  are  v  customers  waiting  in  unit  1.  The  number  of 
arrivals  n  unit  1  between  times  u  and  t  is  The  service 
completion  of  the  customer  in  service  at  t  occurs  at  time  v 
and  there  art  Vg  arrivals  between  t  and  v.  At  time  v  there 
are  v+v^+vg-l  customers  waiting  in  unit  1.  Starting  with  these 
v+Vj+Vg-l  customers  at  time  v,  unit  1  becames  empty  after  r^ 
services  at  time  v^.  Lastly  the  service  completion  of  the 
r+v+v^  customers,  who  arrived  in  unit  1  up  to  time  t,  occurs 
in  unit  2  at  time  Vg.  Now  we  sum  and  integrate  over  all 
choices  of  r,  v,  Vg,  r^,  u,  v,  v^  and  Vg. 

Taking  the  transform  of  (42): 


66 


m 


CO  0»  00 


:*r<s.»>  -  s  r  r  s  sXWV)  „<•) 

“  r^jv^VgOr^N^+Vg  0  2  vi+v2*0 


.  r  .-<**)»  «  r  e-(s+x)v(xv^  m 

J0  V  Jo  V  V  2 


®  »  ®  ®  os  (r)  r+v+v, 

^  Z  TIT  T  gjv  (l)h2  (5)Vv,+v?-l, 

r=0  v*lv1=0VgsCr^=v<-v1+v2-l  1  * 


.  f  e-U+x)t  (xt)  1 


X  dt  j“  .-<*+l>v  ^  dR^v+t) 


Using  lemma  1.1  and  simplifying  the  above  expression  leads  to: 

w 


2XJ  (?’s)  r;3”oo 


"  (r)  r®  -C5+X-Xh2(s)Y  (s)]t 

rsXtufe  2  1  at 

0 

r® -Cs+X-XY..fs)]v  ,  v 

•  I «  1  vr  W> 


Jv  C  X  *> 


dt 


r^Cv=l  *  o 

r®  - us+a-XY., (s)]v 

L  . 

V 


a.  - 


(r)„.  f* 


Sf?(0  J 

r=j 


rv  ■*[?”S+*Y-  ^s)-Xb_(s)Y.  (s)}t  /  \ 

•  J  e  1  1  dt  dHg'(v) 
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.»  -[s+X-XY.  (s)]v 


i 


v  -U-s+XY1(s)-Xh2(s)Y  (s)]t  /  v 

1  2  1  dt  dH«r  (v) 

o 


00 

=  [5-fl+XY1(s)-Xh2(s)Y1(a)]  gjo^^[h2(*+X"XV8^ 


-  hg(£;+X-Xh2(s)Y1(s))J 

+  S  £  8  •  ^  ( ^)h2+v( s ) Y^_1( s )[h  ( s+X-XY  ( a ) ) 
r-0  v=l  Jv  x  L  1  x 

-  h1(l,+X-Xh2(s)Y1(s))J  } 

Using  (3)  and  (1.14)  and  simplifying  further: 

(*5) 

**  r  -i“lf  (s) 

2Xj  (|,8)  =  [|-s+XY1(S)-Xh2(s)Y1(s)J  ^ - ] 

-  qLjl?,h2(«)Y1(s)]  +  [Y1(s)-h1(?+X-Xh2(s)v1(s))] 

•  [(hjjCajYjU))0  -  Yj(§,h2(s))]  ^(s+X-XY^s)) 


-  ^U+X-XhgCs^s))^  | 

This  together  with  vj(?,h2(a))  *  proves  the  lean*. 


Lemma  2,4 


For  R(s)  >  0  and  R(£)  >0  the  transform  2A^  (?>s)  is  given 
by: 

(46) 

2A1  (?,■)  =  [c-s+XY1(s)-Xh2(s)v1(s)J  {^  [«£(?, - - ) 

<  .  _  _1 
-  a"(5,h2(s)\1(s))J  +  [v1(s)-h1(|4X-Xh2(s)Y1(s))] 

•  [b1(s+X-XY1(s))-h1(^X-Xh2(s)Yl(s))] 


[(h2(s)v1(s))i-vi(0 


-  s, 


5+X-s 

X 


where  the  functions  an(-,-)  are  defined  in  (lfi). 
Proof: 

Upon  transformation  of  (40)  yields: 

2Ai  (l,8)  =  oqij)(i)  2Xj  (5»8) 

Substitution  of  lemma  2.3  results  in: 
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(47) 

-~1,  ®  ~  /  \  5_s+^i(s) 

2A"(f,,s)  =  [s-B^CsJ-AhjtslV^s)]  {jJ_0<4  '(5.  J - ) 

•1 

-  0l<”)(S,h2(.)V1(s))}[v1(l.)-h1(^-Ah2(*)Vs»J 

.  [h1(s+X-AV1(«);  h1(5-.X-Xh2(»)Y1(s))] 

•  [(h^Mv^s))1-  e  04|n)(s,o)  ♦  Ji(t,4n,(l,ii2(»)Y1(«)) 

This  and  lemma  2  in  Appendix  B  together  with 

£  4^  (5,0)  =  Yi(^)  prove  the  lemma. 

n=l  ° 

Theorem  2.5 

For  R(s)  >0  and  R(?)  >0  the  transform  ^  (|,s)  of  the 
distribution  function  gWi(t,x)  of  the  total  time  in  the  system  is 
given  by: 

(48) 

=  A*<M)  -  [i  *  *  3^(6..)].  *21. 

where  2A^(S,s)  is  given  by  lemma  2.4. 

Proof: 

Similar  to  the  proof  of  Theorem  2.2. 
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Limiting  Distribution  of  the  Total  Time  in  the  System 

Let  gW(x)  =  lim  gW^(t,x)  and  gO(s)  be  its  L.S.T. 
t  -t  63 

The  existence  of  the  limiting  distribution  can  be  proved  as 
in  Chap.  I  section  5. 

Theorem  2.6 

The  L.S.T.  of  the  limiting  distribution  of  the  total  time 
in  the  system  is  given  by: 

(49)  2u(s)  =  (l-Xc^-XagHl+X  (0,s)]  if  l-XOj-Xotg  >  0, 


=  0  otherwise, 
where  2A^  (o»s)  is  given  by: 


(5°) 

r  -i“^r  00  r  XV  (s)-S 

2^  (0,s)  =  ^-s+XY1(s)-Xh2(s)Y1(s)J  - - - 


) 


-  a^O^U^s))]  + 

[Y^sJ-h^X-Xh^sjY^fl))] 

.  [h1(s+X-XY1(s))-h1(X-Xh2(s}Y1(s))] 

*  [h2(s)Y1(s)-l  +  ^  Can (° ’ h 2 ( S ) Y1 ( 8 ) ) " an ( C’ ’ ~Y  ' )  1  } 

n=l 

and  the  functions  an ( *  s  * )  are  defined  in  (18). 

Proof: 

Equation  (49)  is  obtained  as  in  Theorem  2.3  and  (50)  from 


lemma  2.4  by  taking  £  -*0. 


In  (50)  taking  the  limit  s  -i  0,  ve  observe  that: 


■Jk 


^  f  ^(l-xa^r  +  ax^i-Xc^) 
(l-XaL)2(l-\a1-\a2) 


} 


-  -  {ga-Xa^Cl-t^)2]  (l-Xo^-Xag)}  1 

•[(l-Xa1)(l+Xo!2)(01+02)+2a1a2[(l-XaL)2+X  a2]j 
(by  lemma  1  is  Appendix  B) 

Equation  (52)  together  with  (59)  yields: 


(6o)  ? 

X{(l-XaL)(l+Xa2)(31+B2)+2a1a2[(l-X£zL)  +Xa2]} 

\ 


\ac 


2(1 


6.  The  State  of  the  Ser/er 

In  this  section  we  try  to  answer  questions  of  the  type: 

(i)  What  is  the  probability  that  at  time  t  the  server  is 

busy  (•?  idle)?  (ii)  If  the  server  is  busy  at  t  what  is  the 

probability  that  he  is  serving  in  unit  1  (or  unit  2)? 

(iii)  And  if  he  is  serving  in  unit  1  (or  unit  2)  at  t  what 

is  the  probability  that  he  is  serving  the  r-th  customer  of 

the  cycle  in  unit  1  (*r  unit  2)? 

Finally  we  study  the  limiting  behavior'  the  above 
# 

probabilities. 

r 

We  define  the  following  probabilities: 


1 

j 

1 
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9^r( t/i)  is  the  probability  that  at  t  the  r-th  customer 
of  a  cycle  of  tasks  in  unit  v,  v=l,2,  is  being  served,  given 
that  the  service  started  at  t*0  with  i  >  0  customers  in  unit  1, 
and  its  transform: 


cc 

(61)  9^.U/i)  =  J  e_|t  0  (t/i)  dt,  v=l,2, 

O 

Further,  d^(t/i)  is  the  probability  that  at  t  the  server  is  serving 
in  unit  v,  v=l,2,  given  that  at  t=0  there  were  i  >  0  customers 
in  unit  1.  Then: 


(62)  e  (t/i)  =  e  e  (t/i)  ;  v=i,2, 

r=l  vr 

Lemma  2.5 

The  transforms  evr(^/i)  of  ®vr(t/i)  are  given  by: 


(63) 


I  ,r1J(O[8'r-l)(^l)-/r'1)(?,0)] 


(64) 


>*r(Ul)  -  i  [w.*(4)]  i  i  nr<K,  „ru(0.jB,(fcO)bJ-1(0 

J  =  JL  U=J  — 


*  "T~  ^  oru( ?)€'"'( 5. 0)hr2'l(?) 

where  Qr  (•)  is  defined  in  (1.30),  m^f)  in  (1.52),  g^CV»‘) 
in  (1.6)  and: 
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Ir  =  1  if  r  <  n  fnd  0  if  r  >  n, 
Lr<h] 


Proof: 


Let  r$t!i)  be  the  probability  that  at  time  t  the  server  is  serv¬ 
ing  the  r-th  customer  of  an  arbitrary,  cycle  of  tasks,  in  unit  v, 
given  that  the  queue  has  not  become  empty  in  [O.t]  and  that  at 
t=0  the  service  started  with  i  >0  customers  in  unit  1.  Then: 

(65) 

«.t  /  \  <»t 


oSlr 


(t|i)  -  £  £  £  \  d  Q^n)(u)  J  dG^_l)(u^-u) 

n=0j=l  jx=l  Jo  0  Ju  JJ1 


•  [l  - 


If  the  queue  has  never  become  empty  in  [o,t],  let  the  server  be 
serving  in  the(n+l)th  cycle  of  tasks  at  t,  n  >  0.  At  the  end 
of  the  n-th  cycle  of  tasks  there  are  j  >  1  customers  in  unit  1, 
given  that  at  t=0  there  were  i  customers  in  unit  1,  and  the 
n-th  cycle  of  tasks  ended  between  u  and  u+du.  This  probability 
is  given  by  d  Q^(u).  Starting  with  j  customers  in  unit  1  at 
time  u,  there  are  at  least  (r-1)  services  up  to  time  which 
is  the  last  epoch  of  service  completion  before  t  and  at  the  end 
of  the  (r-l)th  service  >  1  customers  are  waiting  in  unit  1. 
This  probability  is  given  by  the  second  integral.  The  last 
factor  [l-Hj^t-u^)]  ensures  that  at  time  t  the  server  is  serving 
the  r-th  customer.  Finally  we  sum  over  all  choices  of  n,  j, 
u  and  to  obtain  (65). 
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I 


I 


i 


By  a  similar  argument  we  get: 


(66) 


»  eo  as 


oS2r(t  I' '  *  2 


no  j 


■  £  ',i4r'1)<vv[i-vt-2>] 


where  the  factor  I fr<vD  is  preSent  because  if  unit  1  has  v 
services  then  unit  2  can  have  only  at  most  v  services. 

The  Laplace  transforms  of  (65)  and  (66)  lead  to: 


(67) 

,<r(5lD 

(68) 


jf.l  ^!!)  g«l"l)(5) 


aV*i*>  =ih«]  ,=  l  V<v] 


The  usual  renewal  argument  gives: 

(69)  e*r(Mi)  -  o0*rU|i)  +  °1(Ooe*r(t,ll)  ,  v*l,2, . 

Substitution  of  (67)  and  (68)  in  (69)  proves  the  lemma. 
Lenaa  2.6 

If  0^(1  |i)  is  the  Laplace  transform  of  9  (tji),  then: 
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(70)  <(I|D-1  ji0'u(5)[l-^i(6)] 

*  \  “l<5)  £  =V6)  t1  -  vi(?)] 

(71)  e2(?|i)  *  \  *  0rij(e)  [^(0-vJU,h2(5))] 


+  \  .*  or13<5)  h2(S»] 

J— -A. 


where  V  (•)  and  are  defined  in  lemma  1.2. 

Proof : 

This  is  immediate  upon  summing  over  r  >  1  in  (63)  and  (64) 
and  using  lemma  1.3. 

Lemma  2.7 

If  1-Xaj-Xa2  >  0,  the  following  limiting  probabilities 


are; 

(72) 

Vi)  lim  9lr(t|i) 


Xa,  (l-Xa,-Xa2)  I  E  rL1(0)g^"l)(0) 
1  1  6  J-l  1  0  w  J31 


t  -» • 


73) 


(ii)  lim  8,  (t|i) 
t  *t  •* 


(v) 


Xa2(l-XaL-Xa2)  ^  0rlj^gjo 


(0) 


(74) 

(iii)  lim  8  (t|i)  -  Xa  ,  v-1,2, 

t  H  •  V  V 
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Proof: 


The  proof  follows  from  a  Tauberian  Theorem  (Theorem  5> 
Appendix  D): 


(75) 


lim  lim  ie^J%  |i)  ,  v-1,2, 


t  -i  « 


5  -<  o 


vr 


(76)  lim  6  (t|i)  =  lim  (e|i)  ,  v=l,2, 

t  -»®  v  6  -»o  V 
Hence  using  (75)  in  lemma  2.5  we  get  (72)  and  (73). 
Using  (76)  in  (70)  and  (71)  we  have  for  v*l,2: 


(77) 


ter  (1-ter.  -Xa0)  ® 

liB  -  ■  -  T--^— "  i,  •>  or!3(0) 


t  -t  00 


Theorem  4  in  Appendix  B  gives: 


I  J  or  (0)  -  1  ♦  , 
J=1  0  ^ 


(78) 


l-Xo^-Xctg 
1-te^ 


l-Xa^-Xotg 


Substitution  of  (78)  in  (77)  proves  (7*0  • 

Theorem  2.7 

If  0(t|i)  is  the  probability  that  the  server  is  busy 
at  time  t,  given  that  the  service  started  at  t«0  with  i  >0 
customers  in  unit  1,  then  the  Laplace  transform  0#(^|i)  of 
0(t  |i)  is  given  by: 


e*u|1) '  i  j^>  t1  •  'i(5-h2(5!)] 

+  iBl(«^orlj(i)[1-vi[(?,l.2<»)]  , 

Further  the  stationary  probability  that  the  server  is 
busy  is  Xo^  +  Xa 2  and  hence  the  stationary  probability  that 
the  server  is  idle  is  l-Xc^-XOg. 

Proof; 

We  have: 

(80)  e(tli)  -  a^tli)  +  e2(t|i) 

which  gives; 

(81)  e*(!|i)  »  «*(5|l)  ♦  S*(6|i) 

Hence  (79)  follows  from  (8l)  and  lemma  2.6. 

lim  e(tji)  =  lim  e^tji)  +  lim  02(t  |i) 

t  -t®  t-*«  t  H  CO 

=  Xo^  -!•  Xa2  . 

(by  lemma  2.7) 

Expected  length  of  a  cycle  of  tasks.  Starting  with  j  >  0 


customers  in  unit  1,  the  expected  duration  of  a  cycle  of  tasks  is: 
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(82) 


h  q.1(s,l)]8=0  =  “  Is  Vl(s,h2(s))]Ba0 

d(«i+a2) 


l- 


Xa^ 


Expected  sojourn  time.  Consider  the  Semi-Markov  sequence 
{?^,Tn>n  >  0}  defined  in  (1.3).  Let  llj  be  the  expected 
sojourn  time  of  this  process  in  state  j.  Then: 


(83a) 


(83b) 


1  °i  +  a2 

T 7E, 


%  =  , 
J  1 


jfaj+ctg) 


Xa^  ’ 


j  >1  , 


where  (83a)  is  obtained  from  the  fact  that  once  the  process 

reached  the  state  zero,  there  is  a  negative  exponential  idle 

period  with  expected  duration  ^  and  further  a  cycle  of  tasks 

started  with  a  single  cutomer  in  unit  1  whose  expected  value 

is  given  by  (82).  (83b)  is  obvious  from  (82). 

Mean  recurrence  time.  Let  ^  be  the  mean  recurrence  time  of 

J 

state  j  of  the  process  Then: 


lio  -  i  -  Y'(o),  Y( • )  defined  in  (1.1)  , 


(84) 


.  1 

*  xTi-xopXagT 


Let  M.  .(t)  be  the  expected  number  of  visits  to  state  J  by  the 
process  in  (o,t],  given  that  ■  i,  and  m^s)  be  its 


L.S.T.  Then: 
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(85)  *  "i0(s)  xTJ  0^(“>  omiJ<s)-  1  >  1.  3  >  1 

where  Qn  (•)  are  defined  in  (1.30). 

This  is  obtained  from  the  consideration  that  a  visit  to  j 
can  occur  either  with  or  without  an  intermediate  visit  to  the 
state  0  [Neuts  (1969)]. 

From  (85)  we  get  for  j  >  1: 

p,  1  =  lim  s  m  ( s ) 

i  s  -t  0  1J 

3  o^0)  lim  s  mi0(s) 

s  ->  0 

(86)  =  ^(0)  V-'1 

If  j (t )  is  the  probability  that  the  Semi -Markov  process 
is  in  state  j,  given  that  it  started  in  state  i  at  t=0, 
andp  *  lim  p  .(t)  then: 

J  t  h  09 

(QT)  P*  -  ^  ,  i  >  0  , 
i 

Substitution  of  (83),  (84)  and  (86)  in  (87)  yields: 

P*  a  (x  +  tog)  (1  -  ^  -  kxz)  |  (1  -  Xo^) 

and 

p*  -  {X(a1+Q2)(l-XaL-a2)j  o*y(0))  Id-Xo^),  J  >1, 
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7.  Generalizations 

The  Tandem  Queue  With  Mort  Than  Two  Units 

Let  us  consider  m(>  2)  service  units.  The  input  to 
unit  1  is  a  Poisson  process  of  density  X,  and  the  input  to 
the  (r+l)th  unit  is  the  output  from  the  r-th  unit, 
r=l,2, . . .  ,m-l.  After  getting  servi.ce  in  the  m-th  unit  the 
customers  depart  from  the  whole  system. 

At  t=0  a  single  server  starts  serving  in  unit  1.  He 
switches  from  unit  1  either  by  a  zero  switching  rule  or 
by  a  non-zero  switching  rule,  while  he  always  observes  a 
zero  switch  rule  in  all  other  units  2,...,m.  In  all  the 
units  the  customers  are  served  by  the  order  of  their  arrivals 
and  the  server  is  busy  as  long  as  there  is  at  least  one 
customer  in  the  whole  system.  Service  times  are  assumed  to 
be  mutually  independent  positive  random  variables  and 
independent  of  arrival  times. 

Each  cycle  of  tasks  consists  of  m  tasks,  task-1,..., 
task-m,  and  each  busy  period  consists  of  a  random  number  of 
such  cycles . 

Let  H^(« ) , . . . • )  be  the  service  time  distributions  in 
unit  l,...,unit  m  respectively.  The  distributions  of  busy 
period  and  virtual  waiting  time  and  queue  length  are  obtained 
by  replacing  Hg(*)  by  Hg(  •  )*H^(  •  )*. .  • )  in  the  results  of 

two  units.  The  corresponding  moments  are  obtained  by  replacing 

m  m  m-1  m 

a_  by  £  and  09  by  E  0,  +'  2  £  £  a*  a.  * '  * 

i  i-2  1  i»2  i-2  J-i+1  X  J 
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where  av  =  j  xdH  (x)  and  Pv  =  ]  x  dH  (x),  \Fl,...,m 
o  o 

Infinite  Tandem  Queues.  Suppose  that  the  number  of  service 

units  m  is  infinite.  Let  S  be  the  service  time  of  a  customer 

v 

in  unit  v,  v*l,2, .  Then  are  independent  random 

variables  with  distribution  functions  Hg( 

Theorem  2.7.  Convergence  Theorem 
00  00 

(a)  If  a  *  E  a  and  3  =  £  3,.  <  08  then  the  distribution 

v=l  v  v=l  v 

&„.(•)  of  the  service  time  of  a  customer  in  the  first  m  units, 

Hi 

S,  +  ...  4  Sffi,  converges  to  a  probability  distribution  G#(‘) 

with  first  and  second  moments  a  and  3  respectively. 

00 

(b)  The  total  service  time  IS  of  a  customer  converges  in 

v=l 

law  if  and  only  if  for  a  fixed  c  >  o  the  three  series 

(i)  E  f  dH  (x),  (ii)  Ea^  and  (iii)  E  3*^  converge, 
n=l  c  n  n-i  n  n=l  n 

where  =  f  x  dH  (x)  and  3^  =  f  xZ  dH  (x). 
n  J  n  n  j  n 

o  o 

For  the  proof  of  this  theorem  we  refer  to  Feller  (1966). 

If  the  service  time  distributions  are  negative  exponential, 
-u  x 

H^(x)  *  1-e  v  ,  v=l,2,...,  then  by  the  convergence  theorem  the 

distribution  of S,+...+8_  converges  to  a  probability  distribution 
1  m 

*  1 

G  (•)  if  E  -  <».  Q  (t)  gives  the  probability  that  a 
v-1  $ 

customer  will  be  served  in  infinitely  many  units  before  epoch 


t. 
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Under  the  conditions  of  convergence  we  have: 


a?  w 

P  \  E  S  <  x)»  it  *  H  (x) 

Va  v  ’  J  v=i  v 

which  is  the  convolution  of  Hg(. ),.... 

Hence  the  distributions  of  busy  period,  virtual  waiting  time 
and  queue  length  are  obtained  by  replacing  Hg( • )  by 


n  *H  (•)  in  the  results  of  two  units, 
o  v 
v=  2 

Equilibrium  Conditions  of  the  Infinite  Tandem  Queues.  Under  the 
conditions  of  convergence  of  the  total  service  time  of  a 
mstomer,  the  queue  will  attain  its  equilibrium  if 

a> 

1-X  E  a  >  0.  This  follows  from  the  results  of  an  M|o|l 
v=l  v  » 

queue  with  service  time  distribution  it  *  H  (x)  . 

v=l 

The  Tandem  Queue  with  Balking 
Consider  the  tandem  queue  with  two  units.  Let  p  be 
the  probability  that  a  customer  joins  the  queue  in  unit  2 
and  1-p  the  probability  that  he  leaves  the  system  after 
getti.Tg  service  in  unit  1.  The  distribution  of  busy  period, 
virtual  waiting  time  and  queue  length  of  this  model  can  easily 
be  studied  from  the  following  c ons id jrations :  We  assume  that 
all  the  customers  after  getting  service  in  unit  1  go  through 
unit  2  and  get  a  non- zero  service  there  with  probability  p 
and  a  zero  service  with  probability  1-p.  The  distribution  of 
service  time  of  a  customer  entering  unit  2  is: 
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(88)  P  H2(-)  +  (1-p)  U(.) 


Hence  the  distributions  of  busy  period,  virtual  waiting  time,  etc. 
can  be  obtained  from  the  non- balking  case  by  replacing  ( • ) 
by  (88).  To  get  the  moments,  a2  is  replaced  by  p  a2  and  02 
by  p 


In  the  case  of  the  tandem  queue  with  m  service  units, 

let  1?  the  probability  that  a  customer  joins  the  queue 

in  unit  vt]  and  1-p^  the  probability  that  he  leaves  the  system 

after  getting  service  in  unit  v,  v-1, 2, . . . ,m-l.  Here  also 

we  assume  that  each  customer  after  getting  service  in  unit  1 

goes  through  all  the  remaining  (m-l)  units  and  gets  a  non-zero 

v-1 

service  in  unit  v  with  probability  (  s  p . )  and  a  zero  service 


v-1 


i-1 


with  probability  fl  -  (  *  p. )*]  ,  v=2 . m.  The  distribution 

L  i-1  1  J 

of  service  time  of  a  customer  in  unit  v  is: 


(89)  (V  pj  H  (■)  +  fl  -  (V  p  )1  U(-) 

1-1  v  L  i=i  1  J 

Hence  the  different  distributions  of  interest  can  be  obtained 
from  the  non-balking  case  with  m  service  units  by  replacing 
Hv(‘)  oy  (89),  \*-2,...,m. 


8.  Applications 

The  tandem  models  considered  In  Chapters  I  and  II  can 
be  viewed  as  a  modified  alternating  priority  model.  In  the 
alternating  priority  model  [A vi- Itzhak,  Maxwell  and  Miller  (19 65), 
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Neuts  an'".  Yadin  (1968),  Takctca  ( 1968 ) ]  customers  arrive  at  tvo 
service  units,  unit  1  and  unit  2,  in  accordance  with  Poisson 
process  of  densities  X^  and  X g.  A  single  server  attends  to 
two  units  alternately  according  to  zero  switch  rule  and  serves 
the  customers  in  the  order  of  their  arrivals.  In  this 
alternating  priority  model  suppose  that  the  input  to  unit  2 
is  stored  there  as  long  as  the  server  is  serving  in  unit  1. 

Once  the  server  started  serving  in  unit  2  the  input  to  it  is 
shut  off  and  stored  in  unit  1  until  he  switches  back  to  unit  1. 

As  soon  as  the  server  switches  back  to  unit  1  the  stored  input 
of  unit  2  is  released  from  unit  1  to  unit  2.  This  modification 
is  reasonable  in  cases  where  the  arrival  of  a  customer  in  unit  2 
causes  service  interruption  there  or  xn  cases  where  only  those 
customers  of  unit  2  who  have  arrived  during  the  service  time 
of  the  customers  of  unit  1  have  priority  over  the  customers 
arriving  in  unit  1  thereafter. 

The  analysis  of  this  modified  alternating  priority  model 
can  be  easily  deduced  from  our  tandem  model:  Customers  arrive 
at  a  service  system  in  accordance  with  a  Poisson  process  of 
density  X.  Independently  of  others  an  arriving  customer  is  of 
type  1  with  probability  p^  or  of  type  2  with  probability  Pg, 
where  Pi+P2  *  1»  ■  Xp^  Xg  *  Xpg.  All  the  arriving  customers 

pass  through  both  the  units  1  and  2.  A  type  1  customer  receives 
a  non-zero  service  in  unit  1  and  zero  service  in  unit  2,  while  a 
type  2  customer  receives  a  zero  service  in  unit  1  and  a  non-zero 
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service  in  unit  2.  The  service  time  distribution  in  unit  1 
of  an  arriving  customer  is  H^(>)  with  probability  and 
U(‘)  with  probability  p2>  while  his  service  time  distribution 
in  unit  2  is  H2( • )  with  probability  p2  and  U(  • )  with 
probability  p^,  where  U(-)  is  the  unit  distribution.  Hence 
in  our  analysis  in  Chapters  I  and  II  we  replace  H^(x)  by 
PjH^x)  +  p2U(x)  and  H2(x)  by  p^x)  +  p-^x). 


89 


CHAPTER  III 

ALTERNATING  PRIORITY  QUEUES  WITH 
NON-ZERO  SWITCHING 


1.  Concepts  and  Definitions 

This  chapter  discusses  a  queueing  model  in  which  a  single 

server  serves  two  units  1  and  2;  the  input  processes  to  these 

are  independent  Poisson  processes  of  rates  X^  and  X^ 

respectively.  The  server  attends  the  two  units  alternately 

according  to  a  non-zero  switching  rule.  He  continues  to 

serve  in  unit  v  unitl  he  has  given  k  services  without 

v 

interruption  there  or  until  the  queue  becomes  empty  which¬ 
ever  comes  first,  k  v=l,2,  are  positive  integers,  which 

v 

are  called  the  switching  parameters.  The  alternating  priority 
queues  with  zero- switching  have  been  studied  by 

several  authors:  Avi-Itzhak,  Maxwell  and  Miller  (1965), 

Neuts  and  Yadin  (1968),  Takacs  (1968). 

It  is  assumed  that  at  both  units  customers  are  served 
in  the  order  of  their  arrivals .  The  service  times  are 
mutually  independent  positive  random  variables;  independent 
of  the  arrival  times.  Denote  by  H^( • )  and  Hg(*)  the 
distribution  functions  of  service  times  in  units  1  and  2 


respectively. 
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We  use  the  following  notation: 
r°° 

hv(s)  ■  J  e"sxdH  (x),  r(s)  >  0, 

o  v 

oo 

a  =  f  x  dH  (x),  v=l,2. 

V  J  V 

O 

2.  Distribution  of  Busy  Period 

We  recall  that  a  task  is  the  time  interval  spent  without 
interruption  in  a  unit.  A  task  in  unit  v  is  referred  to  as 
a  v-task.  It  consists  of  atmcst  k  consecutive  services 
v=l,2. 

Suppose  that  at  t=0  the  server  starts  serving  in  a  unit. 
The  time  required  for  both  the  units  to  become  empty 
simultaneously  for  the  first  time  is  called  a  busy  period. 

If  the  busy  period  starts  with  the  service  of  a  customer 
in  unit  v  then  the  corresponding  busy  period  is  called  a 
v-busy  period  (or  busy  period  of  type  v),  \j=i,2.  Let 
denote  the  distribution  function  of  type  v  busy  period, 
v=l , 2 . 

The  system  becomes  idle  when  both  the  units  are  empty. 

Hie  idle  period  has  a  negative  exponential  distribution  with 
parameter  +  Xg.  After  an  idle  period  a  new  busy  period 
starts  in  the  unit  in  which  a  customer  arrives  first. 

Remark: 

As  in  Ne^ts  and  Yadin  (1968),  if  the  unit  to  which  the 
server  switches  is  empty  then  we  assume  that  he  instantaneously 


completes  a  task  of  duration  of  zero  there  and  switches  back 
to  the  other  unit. 

Since  the  distribution  of  busy  period  does  not  depend 
on  the  switching  rule,  [We^ch  (1965)],  it  follows  from  Neuts 
and  Yadin  (1968)  that: 

Theorem  3.1 

If  01( s)  and  0?(s)  are  the  L.S.T  of  h^(*)  and  K2(‘) 
respectively  then: 

(i)  For  every  s  wtih  R(s)  >  0,  the  pair  ©1(s)  and  0g(s)  *s 
the  unique  solution  to  the  following  system  of  equations: 

(1)  (a) 

Z1  =  ^(s+^i+^2-^lzl"^2Z2^  ’  Z2  =  h2^S+^l+^2-^lZ''“^2Z2^ 

(2)  (b;  =  ^(s+Xg-XgZg)  >  zg  ”  ^2^®+^l"^l*l^ 

in  the  region  |z.jJ  <  1,  lz2|  <  1,  where  Yy( • )  is  the  L.S.T 
of  the  distribution  of  busy  period  of  an  m|g|1  queue  with 
input  rate  X  and  service  time  distribution  H  (•)>  v=l,2. 

(3) (ii)  01(O+)  =  e2(0+)  =  1  if  and  only  if  l-Xo^-XOg  >  0 
(iii)  If  l-Xo^-XQg  >  0  then  the  means  of  k^*)  and  *2(*) 


are  given  by: 
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3.  The  Basic  Imbedded  Semi -Markov 
Process  and  its  Transition  Probabilities 
We  suppose  that  at  t=0  there  are  >  1  customers  in 
unit  1  and  ig  >  0  customers  in  unit  2.  Furthermore  a  customer 
in  unit  1  is  just  beginning  service.  We  may  also  start  with 
other  initial  conditions. 

Let  us  define  the  sequence  of  random  variables  t^t^,..., 

where  t  =0  and  t  is  the  duration  of  the  n-th  task.  The  odd 
o  n 

numbered  variables  t^,t^,...  are  the  durations  of  tasks  in 
unit  1  and  the  even  numbered  variabxes  tg,t^,...  are  the 
durations  of  tasks  in  unit  2.  Let  ^  be  the 

number  of  customers  in  the  system  (unit  1,  unit  2)  at  the 
end  of  the  n-th  task,  n  >  1  and  =  (i^ig).  Further  let 
b°  a  random  variable  which  akes  values  1  and  2  depending 
on  whether  the  (n+l)th  task  is  a  1-task  *r  &  2-task, 
n  >  1,  CQ=1-  It  then  follows  from  the  regenerative  properties 
of  the  input  and  service  processes  that  the  quadrivariate 
sequence  of  random  variables: 

<5>  K-  il)<  42)-v"2o} 

is  a  Semi-Markov  sequenc°  with  state  space: 

l1’2}  X  l0'1’-"}  x  I0’1’  —  } 

Ta  study  the  transition  probabilities  »f  the  semi-Markov 
sequence  defined  in  (5)  we  define  the  auxiliary*T>Tobability 
functions  and  2o|*^(x)  as: 
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(6a) 


v0«O,(*)  *  5«  U(x)  -  V’1’2- 


and  for  n  >  1,  vG[^(x)  is  the  probability  that,  in  an  m|q|1 
queue  of  input  rate  X^  and  service  time  distribution  H^( • )  the 
initial  busy  period  involves  at  least  n  services,  that  the  n-th 
service  is  completed  before  time  x  and  that  at  the  end  of  the 
n-th  service  there  are  j  customers  waiting,  given  that  there 
were  i  customers  initially,  v-1,2. 

m  r*  *\.y  (\y) 

(6b)  v  ij  (x)  =  J  e  A  Vy)’  v=1’2’ 


(6c) 

r(n+l) 

vGij 


j+i  oc  /  s  -x  y  (x  y)^"r+1 

(x)  =  JvGir  fx-y)e  "  T^iy T-dH  <y),  n>l, 
r=1  •  v=l , 2 , 


Let  ^g^(s)  ^  the  L.S.T.  *f  ^^(x)  and 


(7) 


g[n)(s,t)  =  I  «£>(.)  zj,  |*|  <1,  v-1,2, 
1  v 


j=o 


Then  for  v=l,2: 


(•»*)  =  ** 


(8) 


»  z1'1  hv(s+W} 

U“W>  -  > 
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The  results  analogous  to  lemma  1.1  through  lemma  1.5  are 
easily  seen  to  be  satisfied  by  the  probability  functions 

v9l;’w- 

Let  us  denote: 

(9)  i  -  (ivi2)  ,  (*i»*2) 

Define  the  transition  probability  functions: 

(10)  -  P  {(h-2,!<l)-dx.4!!)-J2.  Tn  ±x 


I  C  =1  E^-i  £^»i  1 

'  Vl  ^Vl  l’Vl  hi 

(U)  Qg(l ,ijx)  >  f{Cn-l.l'1)-J1’?i2)=J2.  \  <  > 


We  have: 

(12) 


h~^2 


Jt  (O  -X9U  (X?u) 

(jVi2i.  • if  ji  2 *•  *x  2 1 


(r)  -V(^)J2'‘2 

•  TJTipr-  ■ if  J1  - 


iz  >  4  2 °> 


Jo  2 


’  U(X)’  lC  ‘l  *  3l  *  °’  l2  -  1 

■  0  ftor  all  other  choices  of  the  indices 

except  for  i^  *  ig  *  0, 
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(13)  Jr1! 

fX  (k2)  ,  (X1U> 

■  3o  az\}2M  •  npyr 


if  d,  >  i»  dx  >  ii  >  °»  i2  > 


Jr1! 


k2  Jt  /,i  -X,u  (!,») 

*  I  d2°i>>  •  -rr-r-rr-  . 


r=i2“°  "  2 


if  dg  =  0,  >  iL  >  0,  i2  >  i 


-  81  j  U(x)  ,  if  i2  =  j2  =  0,  iL  >  1 


=  0  for  all  other  indices  except  ij_  *  *2  “  0 


For  i±  =  0  =  i2: 


jt  -(X1+\g)u 

(14)  ^(O^x)  *  J  ^(liO;  1,  x-u)  e  xi  du  * 

-(\.+X2)u 

(15)  QgfO^x)  3  ]  ft2(°.l>^x-u)  e  X2  du  > 

0 

Let  <i  (i, J , a )  he  the  L.S.T.  of  Q  (i,j,x)  and 

\)  r*S  «"w  V  ^  ~ 

(16)  .  , 

08  “  j-i  J2  , 

q*(i,z,s)  =  Z  E  q  (ij,s)  x  z.  ,  I*  |  <  ,  v=l,2, 

v  ~~  <j1=0d2=ov  1 


From  (12)  to  (15)  we  get: 
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*  i2  r  (ki) 

(i7)  =  *z  {]«<  (S+X2“X2Z2,Z1^ 


*1 


-1 


+  £  ^  (s+X2-X2z2,0)j,  if  ix  >  1> 


r=i1  *1 


-  z22  ,  if  ix=  0,  ig  >  1 
#  f  (^2^ 

<l8>  =  Z1  \2®i  (s+xrxizrz2) 


k2'1 


+  £  ggj  (s+Xj-XjZ^O)},  if  i2  >  1, 


r=i2  2 


-  z^  ,  if  i2  =  0,  i1  >  1 


(!9)  ^(0,z,a)  «  q*(l,0;  z,s) 

1  2 

(2°)  q*(0,z,s)  «  ^2(°>1>  L’s> 

1  2 

Let  R  (i,j,x)  be  the  probability  that  a  busy  period  lasts  for 
n  ~  ~ 

at  least  n  tasks,  that  the  n-th  task  ends  not  later  than  time 
x  and  that  at  the  end  of  the  n-th  task  j  =  (j. ,J„)  customers  are 

1  C 

waiting  i;i  units  1  and  2  respectively,  given  that  the  service 

started  with  i  customers,  i,  >1,  in  unit  1  at  t=0.  Then: 

1  —  ' 


(21a)  *  S.(£»i»x)»  £  x» 
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(21b)  R2n+l^»i»x^  =  E  .  E  r  R2n^~  ~,x“u^ 

\l  =1  v„=0  o 


V1  V 

d<il^i»u)»  n  >  !» 

ix  >1, 


(21c) 


R2n(i/i/x>  “  E  E  f  E2n.l(i/"*x-“) 

V°  V1  0 


dft2^i»u)»  n  >  1>  ix  >  1> 

Further  let  rn(i,j,s)  be  the  L.S.T.  of  R^i^x)  and 

(22)  /(1,1>S)  =  rn(hL,s)  42  , 

X  2 


l*J  <  1.  V  =  1,2, 

h  >1. 


so  that  (21)  gives: 


(23a) 

(23b) 


rl(H/s)  =  , 


00  CO 


W1  A  2~ 


(23c) 


r2n+l(i’i’s)  -  r  £ 

■"I  Vo”  J 

n  >  1,  lx  >  1, 

r2n^~’i>3^  =  E  E 


00  GO 

-  vr=0  ^ 

X  2 


n  >  1,  >  1, 
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(26) 


rn(~~s)  *  2  *l^8>Vi^s)> 

v, =0  v„-l 


and 

(27) 


V1  V° 


4.  The  Queue  length  Process 

Let  us  denote  by  E^(t)  and  ^(t)  the  numbers  of  customers 
who  still  require  some  service  in  units  1  and  2  respectively  at 
time  t.  As  in  Neuts  and  Yadin  (1968)  we  further  denote: 

(28)  =  F1{t1ft)=j1,E2(t).j2 

|et(o).i1,5s.(o)-i2} 

and 

(29) 


where  the  subscripts  1  or  2  demotes  that  at  t  the  server  is  in 
unit  1  or  2  respectively. 

Let  for  v=l,2,  $  (i.j.t)  be  the  probability  that  at  t  there 

V  #s*  As 

are  j  =  (j,  »J-»)  customers  in  units  1  and  2  respectively,  that 
the  queue  is  never  empty  in  (0,t]  and  that  the  original  task 
has  not  ended,  given  that  the  service  started  in  unit  v  at  t*0 


with  i  customers. 
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Further  let  i  he  the  Laplace  transform  of 

V  ^  #v 

$  (i,j»t)  and 

®  ®  J  1  Op 

(30)  '!'  =  E  E  Vi>jl*^*l  z2  » 

Ol=°  j2=0 

K  '  5 


Lemma  3-1 

kh 

The  transforms  ijf  (i,z,0  of  iy  (i,j,t)  are  given  by: 

\)  /V  IV  \J  /W  »•  ■«* 


<1  *  z2  FVV'iVVi1 

-1 

*[*l”hl^+Xl+X2'XlZl"X2Z2^]  } 

*  [l-^  ( l+X^Xg-  X1z1-  X2z2 )  J 

i  +L-  -k  k 

Tl  L1'2!  hl  (5+X1+X2_XiVX2Z2)J 

k  -1  ,  , 

1  f  v-k  le-v  - 

"  E  2l[1-Zl  hl  ^+Xl+X2"XlZl’X2Z2)] 
v=l 

-ig^(?*X2-X2z2,0)}  ,  lj  >1, 

C(i^5)  =  ^{(svVhVVztvV^VVVi-Wl 

(32)  *  ^ g (  ^ ^ ^ 2* Z 1  ” ^  2^ 2'^  J 

■  lz22 


k2"1  _  v-k_  k„-v 


r  2.  2  1 

-  I  z2|  l-Zg  Tig  (^X^+Xg-X^z^-XgZg)  J 

\*=1  w 

.  2g^)(^1-.y1,°)}  ,  i2>i. 
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Proof: 


V,e  have: 


j  - 

V1  jl  f*  (r)  CX,(t-u)l  1 

h(hht)  -  Z  Z  J  V^(u)e  "Tj,'"vTl 


r=0  v=l  o  1 


X2t  ^X2t^ 


^2-i2 


k2_1  “2 


♦2(i;4/t)  -  i  e  j 

r=0  v=l  o  2 


pt  (r\  -X2(t-u)  [X2(t-u)]  2  2 

-far— 


j,-1, 

-X.t  (x,t)  x  , 

*  *  “( j1-i1) •  [1"H2(t-u)J  »  i2  -  lj 


The  probabilistic  arguments  for  these  are  similar  to  those  in 
Chapter  I.  Upon  taking  transforms  in  (33)  and  (3*0  we  obtain: 
(35) 

*1  =  Z22[1'hl^+Xl+X2_XlZl~X2Z2^J^+Xl+X2"XlZl'X2Z2^ 

kl'1 

•  E  {lgir^+X2"X2Z2,Zl^lg!r^+X2'X222,O0’ lH- lj 


r=0  ~  1 


**  i 

*2  =  ^11[l-h2(^X1H2-X1z1-X2z2)](^X1+X2-X1i1-X2z2)- 

V1 

i  {j&i  )(^VV1.z2)  -  2gio^p+xrxiV°)} 


r=0  '  2 


The  lemma  follows  now  by  simplifying  (35)  and  (36)  with  the 
help  of  (l.8b). 


1 

i 

i 

| 


For  r-l,2j  v=l,2,  let  be  the  probability  that 

at  time  t  there  are  j  =  (j., ,j„)  customers  in  units  1  and  2 
respectively,  that  the  queue  is  never  empty  in  (0,t]  and  that 
the  server  is  serving  in  unit  v,  given  that  the  service  started 
in  unit  r  at  t=0  with  i  customers. 


In  terms  of  the  functions  4  and  R  we  have: 

n 


co  ®  oo  »- 

(37)  0^(1^, t)  =  ’^(^.^t)  +  E  E  E  j  '^(v^b-vi) 


v^=l  V2=0  n=l  o 


dR2n(i-^u)>  il^1 


00  CO  CO  pt 


(38)  0-t  9(i  j  j  >b)  =  £  I  I  1  Uv,j,t-u) 

XC  ~  _  r.  >1  C  ~  ~ 


v^=0  Vg=l  n=0  0 


dR2n+l^~~U^  *  il  - 


(39) 


v^=l  Vg=0  n=0  0 


(4°)  022(1^^)  =  + 


®  00  to  -t 


+  E  E  E  j  '■'gC^^t-uJdRgn^^u), 


v^=0  Vg=l  n=l  0 


i2  >1. 
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If  0rv(i>jL>5)  is  the  Laplace  transform  of  (i^t) 


and 


00  00 


‘  ^1  ^2 


(41)  -  f__o  ^ o  »2  . 

1  2 

|*vl  5  r=lj2,  \f1j2, 

then  formulae  (3'7)  to  (40)  give: 

*#  #*  00  00  » 

0*2)  1*^(1,*/ «  -  ♦].  (1,1,5)  +  E  I  Erjj^v,!) 

V1  v2=0  n=l 

K  || 

•  *1  (^0  >  ix  >  1, 


CD  00  00 


(43)  012(i>£>5)  =  s  2  ^  *2n+1^^^2.  > 


v^=0  v2=l  n=0 


il  -  lf 


oo  cm  co 

w  w  w  tj  M 

(44)  *  v£=i  ^ 


'1  2 


i2  >1  , 


(45) 


-¥~y-  y_y 

022^5^)  =  *«>  (**?>$) 


'2  vi>: 


00  OD  CO 


+  *  L,  '2n(^s)+r^5/§)  ’ 


vl=0  V2"'L  n*^ 


i2  -  1  » 

•# 

Denoting  by  n^(i,^t)  the  Laplace  transform  of  n^(i^.t) 
defined  in  (28)  and  (29)  and 
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(46)  «“(i,s,0-  s 


a 3  CD  i  1 

^  V  */4  4  \  jl  J2 

2  L  Vi/A/-'zi  z?  ’ 

d1=o  j2=o  v  ^  1  * 


'z  !  <  1,  v=l,2, 

V  — 


we  obtain: 

Tneorem  3.2 

The  transforms  of  the  joint  distribution  of  queue  lengths 
^(t)  and  ^(t)  and  the  type  of  the  unit  served  at  t  are  given 
by: 

(47) 

■JHfr  x  x  o  ~ - 

*1  (*/£/*)  =  %(*>*/ *)*! U)^' (OtS+^+Xg-X^O^g^ 


+  X2^2i(0A;a/O]i  ijL  >  1> 

and: 

(48) 

**  ■)(•*  ^"1 

4Hfr 

[X1012(1,O;L,O  +  X2022(O,l;z}|)]  ,  ^  >  1 

where  0p^( *,*,*)  are  given  in  (42)  through  (45). 

Proof: 

The  result  follows  from  the  usual  renewal  argument  given 
in  Chapter  I.  For  a  complete  proof  we  refer  to  Neuts  and 
Yandin  (1968). 


105 


5.  Applications 

There  is  a  large  class  of  application  in  which  the  priority 
assignment  follows  more  naturally  from  the  nature  of  the  service 
demanded  than  from  the  urgency  with  which  the  service  is  needed. 
Ir.  many  practical  applications  a  switch  of  service  from  one 
class  of  items  to  another  involves  a  set  up  cost  or  set  up 
time.  The  classification  of  the  input  items  according  to 
similarity  of  service  requirements  is  hence  desirable.  The 
alternating  priority  model  was  first  discussed  by  Avi-Itzhak, 
Maxwell  and  MUle"  (1965).  They  considered  the  alternating 
priority  model  with  zero  switching.  The  model  wa  considered 
in  Chapter  III  is  the  non-zero  switching  case  which  is  a 
generalization  to  zero- switch.  Although  the  analysis  of  the 
non-zero  switching  model  is  very  complicated,  it  is  more 
practical.  In  the  case  of  a  device  controlling  traffic  at  an 
intersection  the  zero  switch  rule  allows  one  stream  of  vehicles 
access  to  the  intersection  as  long  as  there  are  vehicles  in 
this  stream  and  a  steady  input  of  vehicles  in  this  stream 
delays  other  streams  indefinitely.  A  compromise  rule  is  to 
allow  a  certain  number  k^  of  vehicles  of  one  stream  access 
t«  the  intersection  and  then  that  stream  is  stopped  and 
to  allow  a  certain  number  k^  from  another  stream,  etc.  The 
optimum  numbers  and  k^  may  then  depend  on  traffic  conditions. 
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CHAPTER  IV 

A  PRIORITY  RULE  BASED  ON  THE 
RANKING  OF  THE  SERVICE  TIMES  FOR 
THE  M|G|l  QUEUE 

1,  Concept?  and  Definitions 
This  chapter  presents  mainly  the  content  of  the  article 
by  Nair  and  Neuts  (1969).  Here  we  propose  a  priority  rule 
based  on  the  length  of  service  demanded  by  a  customer. 

Takacs  (1964)  discussed  a  priority  queue  based  on  the  rankings 
of  the  service  times  of  the  customers  and  obtained  the 
asymptotic  moments  of  the  virtual  waiting  time  assuming  that 
a  customer  with  a  shorter  service  time  has  priority  over  a 
customer  with  longer  service  time.  Here  we  consider  a 
different,  but  related  problem. 

We  first  recall  a  branching  process  description  of  the 
M|G|l  queue  suggested  by  Kendall  (1951)  and  investigated  by 
Neuts  (1969).  Suppose  that  at  time  t=0  there  are  i  >  1 
customers  in  the  queue  and  that  one  of  them  is  just  entering 
service  at  that  time.  These  customers  form  the  first 
generation  and  their  total  service  time  is  the  lifetime 


of  the  first  generation.  Customers  arriving  during  the 
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lifetime  of  the  first  generation,  if  any,  make  up  the  second 
generation,  with  its  lifetime,  and  so  on.  If  at  the  end  of 
the  first  or  a  subsequent  generation's  lifetime  there  are  no 
customers  in  the  queve,  then  there  is  an  idle  period  at  the 
end  of  which  a  customer  arrives  who  makes  up  the  first  generation 
of  a  busy  period. 

It  is  clear  that  the  life  time  of  a  generation  doer  not 
depend  on  the  order  in  which  customers  have  been  served  dur¬ 
ing  it .  We  will  study  the  virtual  waiting  time  for  the  M  |g  |l 
queue  under  the  assumption  that  within  each  generation  customers 
are  served  in  the  order  of  shortest  (or  longest)  service  times. 

We  will  call  these  policies  the  shortest  processing  time  (SPT) 
and  the  longest  processing  time  (LPT)  disciplines,  respectively, 
and  compare  them  to  the  first-come,  first-served  (FCFS) 
discipline.  Once  the  rearrangement  is  achieved  within  a 
generation,  the  incoming  customers  thereafter  do  not  upset 
it;  hence  the  question  of  service  preemption  does  not  arise 
here. 


2.  The  Basic  Imbedded  Semi-Markov  Process 

We  assume  that  at  t=0  there  are  i  >  0  customers  in  the 

queue  and  that  the  one  with  shortest  service  time  enters 

service  immediately.  A  sequence  of  random  variables 

T  is  defined  as  follows:  T  =0  and  T  is  the  time 

o  1  c  n 

at  which  all  customers,  if  any,  present  at  Tr  ^  complete 

service;  if  there  are  no  customers  at  T  . .  then  T  is  the 
’  n-1  n 
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time  which  the  first  customer  who  arrives  after  T  . 

n-1 

completes  service.  That  is,  Tn  is  the  time  of  service  completion 
of  n-th  generation,  if  the  n-th  generation  is  not  empty.  On  the 
other  hand,  if  the  n-ah  generation  is  empty,  then  Tr  is  the  time 
of  service  completion  of  the  first  customer  who  initiates  the 
first  busy  period  after  time  T 

Let  ^(t)  denote  the  number  of  customers  in  the  system, 
at  time  t+0,  who  still  require  some  service.  Then  the 
bivariate  sequence  of  random  variables: 

(1)  ltd,).  Vi  ■  V  “  2  °1 

is  a  Semi-Markov  sequence. 

We  define  the  taboo  probabilities: 

(2)  U(x), 

and 

oQij)(x)  =  ^n-*’  5(Tn)=j,|(Tv^0,v=l,2,...,n-l 

1  5(T0)*i]  ,  n  >  1 

3.  The  Virtual  Waiting  time  Process 
Consider  an  M|G|l  queue  that  has  a  Poisson  input  with 
parameter  X  and  a  continuous  service -time  distribution  function 
H(*)  with  finite  mean  a.  We  denote  by  Tj(t,x)  the  waiting  time 
of  a  virtual  customer  arriving  at  t  whose  service  time  is 
x  >  0,  where  the  M  |G  |l  queue  observes  an  SFT  discipline,  and 
T(t,x)  the  corresponding  virtual  waiting  time  in  an  m|g|1 


109 


queue  with  an  LPT  discipline.  Let 

(3)  Wi(t,x,y)  =  ?\o  <  Tj(t ,x)  <  y  |  5(0)  =  i}  , 

(4)  Ai(t,x,y)  =  p|o  <  J}(t  >x)  <  y,  T}(t,x)  j  0 

for  all  i€(0,t]  |  5(0)  =  i} 

Then  as  in  (l.6l)  we  have: 

ft 

(5)  W^t^y)  =  Ai(t,x,y)  +  J  A^t-x^y)  dM^-r) 

A  o 

+  p{3(t,x)  =  o  I  5(0)  =  i}u(y) 

Let  Wi(t,x,s)  and  Ai(t,x,s)  respectively  be  the  L.S.T.  of 

W.JtjXjy)  and  A^tjXjy)  with  respect  to  the  variable  y  and 
**.  .  **.  . 

let  W.  (5jX,s)  and  (5>x,s)  respectively  be  the  Laplace 

-x* 

transforms  of  Wi(t,x,s)  and  Ai(t,x,s)  with  respect  to  t. 
Further  we  denote: 

(6)  H(s)  =  if  0  <  z  <x  , 

=  0  otherwise, 

and  h(s),  h(s),  Qq^(s)  the  L.S.T.  of  H(-),  H(0  and  0Q^(-) 
respectively,  and 

(7)  q.[n)(s,2:)  =  Z  ,<”>(.)  zj  ,  |«|  <1. 

j=o 


no 


Lemma  4,.l 

## 

Po:r  R(f)  >  0  and  R(£)  >  0,  the  transform  (^,x>s)  of 
Ai(t,x,y)  is  given  by: 

(8)  Aj*(S,x,s)  =  ^  [h^(§,Z)  -  h^^z')]  > 

where : 

(9)  Z  *  h[s+XH(x)[l-h(s)]},  Z'  =  h[f>+XH(x)[l-h(s)]}  , 


and 

(10)  bo(5,z)  =  z,  hn(5,z)  =  h[^+x-xhn_1(5,z)]  ,  n  >1  , 


Proof: 


We  have: 


(11) 

A1(t,x,y) 


(u)(v) 


n=0  j=l  v*0 


do*ij 


(n)(u)e'X(v"u) 


£  (?)  Hk(x)[l-H(x)f“1C  H*k)(y-v+t)dvH(j'(v-u)  , 
k=0 

where  R(*)  is  defined  in  (6),  H^(*)  and  H^(')  are  the  in¬ 
fold  convolutions  of  H(*)  and  H(-)«  The  probabilistic 
argument  to  get  (11)  is  the  following:  If  the  queue  has 
never  become  empty  in  (0,t]j  let  the  last  beginning  of  the  life 
of  a  generation  occur  between  u  and  u+du  and  let  there  be  j 
individuals  in  that  generation.  Let  the  end  of  the  life  time 
of  that  generation  be  between  v  and  v+dv  (v  >  t).  In  the 
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internal  (u,v),  v  >  0  customers  arrive,  and  they  have  priority 
over  the  virtual  customer  if  and  only  if  their  service  time 
does  not  exceed  x.  If  there  are  k  such  customers,  0  <  k  <  v, 
then  the  distribution  of  their  total  service  time  is 
The  formula  (11)  is  obtained  by  using  the  independence 
properties  L-id  summing  ever  all  allowable  values  of  n,  j, 
v,  k,  u  and  v. 

Taking  the  transforms  of  (ll): 


(12) 


Aj,(t,x,s)  = 


t  00 

'8tJ  I 

o  Jt 
(u)(v) 


n=0  2=1 


dy  H(j)(v-u) 


(13) 


\  u,x,s) 


1 


.^(S)  \  •5Co+XH(x)(l-STCa))] 


n=0  j=l 


-  hJ[tnH(x)U-h(s))]} 


(«,*')} 


9 


where  Z  and  z'  are  defined  in  (9). 

Now  the  lemma  follows  from  lemma  1  in  Append!-  C. 
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Theorem  4.1 

For  R(s)  >  0  and  R(  0  >0  the  transform  W  (fc,,x,8) 
of  the  distribution  function  of  the  virtual  waiting  tiue 
Tj(t,x)  is  given  by: 

W  \  (e,,*, *)  =  -^Ty  £  [^(5,2)  -  h^z')] 

*  Fx^lrl1  *  tfit  ) 

wht*e  Z  and  Z;  are  given  by  (9)  and  h  (•,.)  by  (10). 

Proof: 

Taking  the  transform  of  (5/  »e  get,  as  in  Theorem  1.2, 

that: 

(15)  Wt  (§,x,s)  =  (5,x,a;  +  g+xlffilj  |_1+XA1  (S»X»8)J 

The  theorem  now  f.llows  frem  lemma  4.1. 

Limitiag  Behavior  o*  Virtual  Waiting  time  Process 

Let  W(x,y)  =  lim  W. (t,x,y).  The  existence  of  this 
t  h  <*> 

miting  distribution  can  be  demonstrated  as  in  Theorem  1.3* 
Theorem  4.2 

The  L.S.T.  u(x.s)  of  the  limiting  diatiibution_W(x,y)  of 
the  virtual  waiting  time  T)(t,?)  is  given  by: 

to  _ 

(16)  <a(x,s)  -  (l~Xa){l  -  £  Z  [hn(0,z)-hn(0,z) j  },  if  l-*a  >  0, 


0  otherwise, 
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where  hn(*>0  are  given  by  (10)  and 

(17)  Z  =  hjs+XH(x)£l-Jr(s)Jj-  ,  %  =  h  |aH(x)  jjL-?T(t>)J  j 
Proof; 

Similar  to  the  proof  of  Theorem  1.3.  As  in  (l. 65)  we 
obtain: 

(18)  cu(x,s)  =  (l-Xcr)  [1  +  X  (0,x,s)]  if  1-X  a  >  0 

=  0  otherwise, 

Substitution  of  lemma  4.1  in  (18)  proves  the  theorem. 

Taking  the  limit  as  s  0+  in  (l6)  we  observe  that 
oj(x,0+)  =  1. 

The  Moments  of  the  Limiting  Distribution 
We  use  the  following  notation: 

00 

0  “  I  u  dH(u)  , 

Jo 

00 

V  =  J  u3  dH(u), 
o 

rx 

a  =  u  dH(u)  , 


Vx  *  J*  u3  dH(u)  , 


and 

*n(x<S)  =  hn(0,Z)  -  hn(0,Z)  , 


(19} 
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where  Z  and  Z  are  given  in  (17 ) . 

In  terms  of  the  functions  *  formula  (16)  yields: 

00  \Jf  (x  5  s ) 

(20)  tu(x ) s )  -  (l-xa)  {1  -  x  I  -  ]  , 

L  n=0  9  J 


By  lemma  3  in  Appendix  C  the  series  E  ^(xjs)  is  dominated 

n=l 

by  a  convergent  series  if  1-X  a  >  0.  Hence  by  Lebesgue 
dominated  convergence  theorem,  term  by  term  differentiation 


gives: 


an  “  <(x,s)  -  *  (x,s) 

rr  «(x,8)l  =  Xd-to)  E  -2 - , - 2 - .1 

38  4=0  n=0  S2  4=0 


Applying  l'Hopital's  rule  twice  we  get: 


S  *!(x,o) 
n=0 


where  the  number  of  primes  denotes  the  number  of  derivatives 
taken  in  succession  with  respect  to  s.  Similarly: 

(23)  . .  Mil isl  l  *>,0) 

ds^  Js=0  3  n=0  n 

From  (19)  we  have  for  n  >  0, 


<(-■0,  ■  g)  l  *  y  i  ,=0 


NlA-*'" 
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(27)  E  <l(x,0)  = 

n=0  n 


-1 


(l-Xa)(l-\2a2) 


+  (l-xV)'1 


(l+3Xax+3X2a2)[Y(l-\2a2)+3X2a32]} 


Let  M^(x)  and  V^(x)  denote  respectively  the  first  and  second 
moments  of  the  limiting  distribution  of  T)(t,x).  Substitution 
of  (26)  and  (27)  into  (22)  and  (23)  respectively  leads  to: 

X  0(1  +  2X  a) 


(28) 


and 


v 


x)  = 


2(1  -  X2  a2) 


(29)  v  (x)  =  - — -  {3X0B  +  (l-X3a3)“1(l+3Xav  +  3X2a2) 

$  3(l-xV)  L  X  xx 


.  [v(l-X2a2)+3X2a  32]  } 


4.  The  Longest  Processing  Time  Discipline 
In  the  longest  processing  time  (LPT)  discipline,  within 
each  generation  the  customers  are  ordered  according  to  their 
length  of  service  times,  with  highest  priority  going  to  the 
customer  with  longest  service  time.  The  virtual  waiting  time 
process  of  the  present  case  can  be  treated  as  in  the  case  of  SFT 
discipline.  As  we  have  denoted  ‘n(t,x)  is  the  virtual  waiting 
time  of  a  customer  arriving  at  t  whose  service  time  is  x  >  0 
in  the  case  of  LFT  discipline.  The  Laplace-Stielt jes  transform 
of  the  limiting  distribution  of  T,(t,x)  can  be  obtained  as  in 
(16): 
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(30) 
where 

(31) 


and 

(32) 


(l-Xa)jl  -  |  £  [hn(0,0  -  hn(0,C)]},  for  1-Xa  >  Os 


C  =  h(s  +  X[l-H(x)]  [1  -  h(s)]}  , 

C  •  h{X[l  -  H(x)]  [1  -  h(s)]}  , 

•«  _ 
h(*)  is  the  Laplace-Stieltjes  transform  of  H( • ) 


rs.  . 
H(z) 


H(z)  -  H(x) 
1  -  h(x) 


if  z  >  x 


=  0  otherwise  , 

The  first  and  second  moments  of  the  limiting  distribution  of 
T](t,x)  are  obtained  from  (30)  as: 


XP(1  +  2  X  a) 

(33) 

2(1  -  X  a  ) 


and 


(34)  v— (x)  =  Vy- 
11  3(1-XV) 


^3XPB*  +  (l-xV)'V+3Xa*+3X2a*  ) 

[v(l-X2a2)  +  3xW]  ) 


where 


•if 


a 


X 


a  -  a 

x 


5.  Comparison  of  the  SFT,  LPT  and  FCF5  Disciplines 
Let  T|(t)  be  the  virtual  waiting  time  of  a  customer  arriving 
at  t  in  an  m|G|1  queue  with  FCFS  discipline  and  let  be  the 
mean  of  the  limiting  distribution  of  T|(t).  Then  it  is  Known  that 
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^  *Si  2(l-Xa)  * 


From  (28),  (33)  and  (35)  we  observe  an  interesting  relation¬ 


ship  among  M^(x),  M^(x)  and 
(36)  =  \  [m3(x)  4  m-(x)] 


Also, 


(37)  Mn(x)  <  <  M^(x)  if  and  only  if  <  |  , 

Again,  M^(X)  and  M^(x)  are  random  variables  with  respect  to  X, 
which  has  a  distribution  function  H(0*  If  we  denote  by  Ex 
the  expectation  with  respect  to  the  random  variable  X,  then: 
r00 

Ex  Mg(X)  =  J  M^(x)  dH(x) 


- T1  +  2ta~  [  u  H(u)  dH(u)1  , 

2(l-X  a  )  L  Jo  J 


CO 

(39)  Ex  M^X)  =  [l  +  2X  J  u  H(u)  dH(u)]  , 

In  particular  if  H(x)  =  1  -  e  ,  x  >  0,  then: 


HgOO 


2«1-P) 


E*  *vx) = 
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where  p  is  the  traffic  intensity  ^  .  Hence  in  the  case  of  an 
M|m|i  queue: 


(43)  Ex  Mjj(X)  <  <  Ex  M^(X) 

Further  examination  reveals  that: 

A  +  2\  a  v 

w  Vx) '  (xrnr )  Mn 


That  is,  the  steady  state  expected  virtual  waiting  time  under 
SPT  rule  is  obtained  from  the  steady  state  expected  virtual 
waiting  time  under  FCFS  rule  by  multiplying  with  the  factor 
(1  +  2X  a  ) j  (1+  kx)  which  increases  monotonically  from 
l/(l  +Xa)  for  x  -  0  to  (1  +2Xa)/(  1  +  \a t)  for  x  ».  The 
factor  (1  +  2Xa  )/(l  +  Xa)  >  l/(l  +  Xa)  >  \  ,  since  Xa  <  1 
by  the  steady  state  condition.  Hence: 


for  all  x  >  0, 


and  M^(x)  1  \  as  ^  ^  aod  for  small  x.  Again, 

(l+2Xax)/(l+Xa)  <  (^~)  =  ?  "  l/(l+Xa)  <  |  which  implies 


that : 


M^(x)  <  |  for  all  x  >  0, 


and  M^(x)  t  |  as  Xa  t  1 


ana  for  large  x  Thus; 


C*5)  |  5  «j(x)  <  \ 


Similarly  we  observe: 
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1  *  2  \  a 

(46)  MTitx)  ■  ■ -~a  X-)  S  • 

If  we  draw  the  graphs  of  y  =  M^(x)  and  y  =  M^(x),  then  it 
is  easily  seen  that  they  are  symmetrically  situated  on  either 


side  of  the  line  y  = 


Hence  whenever 


M^(x) 


satisfies 


1 - ! - - - 1 - ? 

o  Ay 

Figure  3. 

The  Comparison  Graph 

the  inequality  (45)  N^(x)  also  satisfies  the  same  inequality 
but  realizes  in  the  reverse  direction.  They  are  concurrent 
with  when  ax  =  This  ir  graphically  shown  in  Figure  3- 
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6.  A  Renewal  Argument 
Equation  (11)  can  be  written  as: 


(47) 


t  °°  t-u+y  / .  \  00  /  \ 

Ut,x,y)=  f  £  R  (du)[  HU'(dz)  £  e’Xz 


'0j  =  l 


v=0 


k=0  K 


which  is  obtained  by  replacing  v-u  by  z  in  (ll)  and  defining 

03  /  \ 

£  Q:n'(t)  =  R.  .(t).  R. .(t)  is  the  expected  number  of  times 

n=o  0  13  o  ij  o  xj' 

state  j  is  entered  without  visiting  the  state  zero  in  [0,t], 
sbarting  at  state  i.  Defining: 

(48) 


F,(t-u,x,y)  .  f W>  S  ^  ^  l  „*W 
3  Jt-u  v=0  ■  k=0  K 

•  [l  -  H(x)f  H^k;(t+y-u-z)  , 


we  rewrite  (47)  as: 

t  « 

(49)  A  (t,x,y)  =  E  R  (du)  F  (t-u,x,y)  , 

1  o  j,l°  iJ  J 

(49)  together  with  (5)  gives: 

ft  « 

(50)  W.(t,x,y)  =  j  £  qR .  (du)Fj(t-u,x,y) 

0  j~l 

rt 

+  J  ^(t-r,x,y)  dM^x) 

o 


+  P(3(t,x)  -  0  |  $(0)  =  i]  U(y) 


By  Smith's  Key  Renewal  Theorem  {Theorem  k,  Appendix  D)  and 
lemma  1.8: 

(51)  W(x,y)  =  lim  W  (t.,x,y) 

t  -S  <o 

<0  CD  03 

=  |  £  ~  F,(T,x,y)dT  +  M  A.  (t ,x,y)dr 

Jo  j=l  J  »  Jo  1 

J 

+  (l-Xa)  U(y)  ,  if  1-X  a  >  o  , 

=  0  otherwise 


where  p  is  the  mean  recurrence  time  of  state  j  without 

°i 

visiting  state  0  and  p  is  the  mean  renewal  time  of  the  general 
renewal  process  formed  by  the  beginning  of  busy  periods. 

From  (1.52)  and  (2.86): 

p  1  -  0  and  p  1  =  X(l-Xa) 

°j 

Hence  further  simplification  of  (51)  gives: 

(52) 


GO  CO 

W(x,y)  =  (l-Xa){u(y)+X  Z  R..(+®)[  F .(x,x,y)  dx  1 
L  j=l°  iJ  o  J  J 


if  l-xa  >  C, 


=  0  otherwise, 
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7*  An  Exact  Comparison  of  the  Waiting  times 
Under  Three  Priority  Rules 

A  number  of  compariso/.s  between  SPT,  LFT  and  FCFS  rules 
were  carried  out  in  the  earlier  sections  in  regards  to  expected 
waiting  times  in  the  equilibrium  state.  Many  questions  which 
involve  more  than  expected  values  may  be  asked  however  and  in 
order  to  answer  them  an  exact  comparison  of  the  waiting  times 
as  random  variables  needs  to  be  made. 

We  may  "visualize"  the  definition  of  the  three  random 
variables  T]( t , x ) ,  T|(t,x),  T](t)  on  a  common  probability  space 
as  follows.  Imagine  that  a  customer  joining  the  queue  at  time 
t  consists  of  three  identical  parts  1,  2,  3  all  requiring  a 
processing  time  x  >  0.  Part  I  waits  in  front  of  a  server  operating 
under  the  SPT  rule,  part  2  in  front  of  a  server  operating  under 
the  LPT  rule  and  finally  part  3  waitr  in  ftont  of  a  unit  governed 
by  the  FCFS  rule.  Then  Tj(t,x),  Tl(+  )  and  'H(t)  are  the  waiting 
times  of  parts  1,  2  and  3  respect  2ly. 

Ar.  Auxiliary  Calculation 

Consider  the  time  points  t  and  t+t ' ,  t  >  0,  t '  >  0.  The 
probability  that  during  the  interval  t,  customers  arrive 
whose  service  time  is  less  than  x,  j^  whose  service  time  is 
greater  than  x  und  that  during  (t,t+t*)  j^  anu  arri 'f.  with 
service  times  respectively  less  and  greater  than  x  is  given  by: 
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We  assume  that  x  is  a  ooint  of  continuity  of  H( < )  so  that 
the  probability  that  one  or  more  customers  have  service  time 
exactly  equal  to  x  is  zero. 

Let  II*  and  I)g  be  the  total  service  time  of  all  customers 
in  (0,t)  with  service  time  respectively  less  and  greater  than  x. 
Similarly  and  are  the  corresponding  quantities  for  the 
customers  arriving  in  (t,t+t/). 

We  define  W(t,t*  ;  x^x^x^x^)  as  the  probability  that 
for  given  t  >  0  and  t '  >  0,  the  random  variables  U^,  U*,  U*, 
satisfy: 

Ug  <  Xg,  U*  <  x3>  <  x^  . 


It  follows  from  (53)  that: 

(54)  W(t,t';x1>x,,x  ,x.)  =  E  E  E  E  e 

Jl=°  j2=0  j.3=°  V" 

i,  ,  J,  .  i 

"  I  i 

'3 


®  00  00  CD 


-xt-xt ' 


I»  H(x)]  1  p.t'u(x)]  l  tXtCl-H^)])' 

•*1*  ^3*  ^2" 

■  HUl\x1d>3\^h\xr)^h' 


M 
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where  !T  and  IF  are  defined  in  (6)  and  (32)  respectively.  Upon 
taking  Laplace-Stielt jes  transforms: 


(55)  W*(t,t/;slis2Js3Js^) 


~8lxl“s2x2“s3x3-s4x4 


we  obtain: 


>x4 


W(t,t  >xi’^2,X3,x4^ 


(56)  W  (t,t/;«1,s2>svs.)  -  2  £  2  2  e 

Jl“°  V0  j3=0  V0 


=  expj-Xt-Xt  VxtH(x)?T(s^)+XtCL-H(x)35l,(s2) 

+  Xt /H(x)h(s3)+Xt  ,[1-H(x)]h(s4)) 

We  now  return  tc  an  K |G jl  queue,  which  we  consider  at  time 
t.  We  define  the  following  five  random  variables.  U  is 
the  length  of  time  beyond  t  until  the  generatior  of  customers 
in  service  at  time  t  completes  its  service.  and  U2  axe 
respectively  the  total  service  times  of  the  customers  with 
processing  times  less  than  and  greater  than  x  who  have  joined 
the  queue  since  the  beginning  of  the  service  time  of  the  current 
generation  and  before  t.  U7  and  U4  are  respectively  the  total 


service  times  o Z  the  customers  with  processing  times  less 
than  and  greater  than  x,  who  join  the  (lueue  during  the  time 
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interval  (t,t+UQ). 

If  at  time  t  the  server  is  idle  all  the  five  variables 
are  zero. 

We  shall  express  the  joint  distribution  of  the  waiting 
times  TJft,x),  ^(t,x)  and  TJ(t )  in  terms  of  the  joint  distribution 
of  the  random  variables  U^,  j«0,...,4. 

The  Joint  Distribution  of  U. ,  j=0,...,4 

Let  oR^(t,xo,x^,x2,X2,x^)  be  the  probability  that  in 

(0,t)  the  queue  has  never  become  empty  and  that  the  variable 

Uj,  j=0,...,4  associated  with  the  time  point  t  satisfy 

U  <Xj,  j=0,,..,4,  given  that  at  t=0  there  were  i  >  1  CwStomerc 

in  the  queue,  one  of  who  was  beginning  his  service  at  that  time. 

Then;  x 

(57)  oRi(t,xo,x1,x2,x  ,x4)  =  S  I  j  J  do^)(x) 

n=0  v=l  o  o 
(t#)  (t) 

•  dH^(t+t'-r) 

♦  W(t-T,t';x1,x2,x3>x4) 

The  probabilistic  argument  for  this  is  the  following:  At 
some  time  ?  prior  to  t,  the  generation  in  service  at  t  enters 
service.  There  are  some  number  v  >  1  customers  In  it,  so 
that  the  duration  of  the  total  service  time  of  these  v  customers 
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has  as  its  distribution  the  v-fold  convolution  H^v\*)  of  H(*)* 
If  U  <  xq  must  hold,  then  the  total  service  time  of  these  v 
customers  cannot  exceed  t+xQ.  The  other  requirements  <  x, , 
^2  -  x2’  U3  -  x3i  x4  account  for  the  factor  W(t-r; 
t»xi»x2>x3»x4)*  ®e  probabilities  0Q^(*)  are  defined  by  (2). 
Taking  L.S.T.  of  (57): 


.»  «“  ,»  .* 


ah(«.vws)-I  •-**!  !!il 


O  0  o  o  o 


-S^-S^-SgXg-S^X^-S^ 


dx0,x1,x2,X2»x^  oRi^t,xo'xl,x?,x3,x4^ 

and  recalling  (56),  we  obtain: 

(59) 

A  (5*S0>S1,S2,S3’V  =  Z  2  q^(|)j  J  exp^-§t^-Xt^ 

n=Ov=i  o  o 

(V'V 

"  *V  SoXontlH(x^  S1 ) ' +*t1  [l-H  (x )  ]fT(  s  2 ) 

+  XxoH(x)h(s;J)>Xxo[l-H(x)]h(s4)JdH^^(t1+xo)dt1 
=  {^so-XH(x)D^81-HIs3)]-X[l-H(x)]EF:(s2)-l?ts4)]}“1 

S  1  l[n)  k,h[^«0-kH(x)h(.,)-X(l-K(x))et.J]) 

-  oq[n)U,hU+X-XH(x)fr(5l)-\(l-H(x))h(s2)]}  j 
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in  terms  of  the  functions  oq|n^(sj,z)  defined  in  (7). 

Next,  let  R^tjX^x^x^XyX^)  be  the  probability  that 
at  time  t,  the  random  variables  associated  with  t  satisfy 
Uj  <xj>  j=0,l,2,3,4,  given  that  at  t=0  there  were  i  customers 
in  the  queue. 

The  standard  regeneration  argument  as  in  (l. 3°)  leads  to: 

(60)  R^tjx^x^xgjx^x^)  =  oRi(t,xo,x1,x2,x3,x1|) 

+  J  oH1(t-u,::0,x1,x2,x  x4)  dM^u) 
o 

+  ?{t(t)=0  |  5(0)=i]  U(x1,x2.v3,x4) 

where : 

U(x1,x2,x3,x4)  =1  if  x^  >0  for  all  j=l,2,3,4 
=  0  otherwise. 

Upon  taking  transforms  in  (60)  we  obtain  as  in  (15): 

^  H  ||  || 

(61)  R^  (  Sq  > Sg , S^jSj^)  -  Sp, 3^ , Sg, Sj^) 

+  y1  ( t) C§+x- xy( %) ]_1  {i+x  A"(^0.W3.S)} 

When  1-Xa  >  0,  the  existence  of  a  joint  limiting  distribution 
for  Uj,  j«0,l,2,3,4  is  guaranteed  as  in  Theorem  1.3.  Further 
when  1-Xor  <  0,  R^tjX^x^Xg.x^x^)  tends  to  zero  for  all  i  and 
xj  Z  j=0,...,*».  Since  the  limiting  distribution  exists  when 
1-Xot  >  0,  its  transform  is  given  by: 
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(62) 


s1,s2,s3>su)  =  ^lim^S  Ri  (^s0>a1»fi2’83}S4^ 

f  "i 

=  (1-ta)  ]l+\  oP1  (0+,8o,81,S2,83,8u)j 


The  Joint  Distribution  of  Tj(t,x),  T](t)  and  T](t,x) 

The  random  variables  Tj(t,x),  ^l(t,x)  and  1|(t)  are,  for  each 
t  >  0,  related  to  the  random  variables  U^,  j=0, 1,2,3, 4 
associated  with  the  time  instant  t  by: 

(63)  j)(t,x)  =  U,  +  Ux  +  U3 

Tl(t)  =  Uo  +  Ux  +  U2 

Ti(t,x)  =  UQ  +  U2  +  Uu 

That  till 8  is  indeed  so,  we  argue  for ,T](t,x).  The  other  cases 
are  similar.  Consider  a  virtual  customer  with  service  time  x 
arriving  at  time  t.  He  has  to  wait  until  all  customers  of  the 
present  generation,  if  any,  have  been  served.  This  is  a  length 
of  time  Uo<  Next,  in  the  next  generation,  all  customers  with 
service  time  less  than  x  are  served  ahead  of  him.  Regardless 
of  the  actual  order  of  service  the  total  amount  of  processing 
time  required  by  all  customers  with  service  time  less  than  x 
is  +  U3-  is  the  processing  time  of  those  who  preceded 

him  and  U3  that  of  those  who  succeeded  him  in  the  arrival 
sequence.  He  have: 
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4Hf 

=  Ri  (5,C1+C:+C3,Cl+C2>C2H3,C1,C3) 

where  (•,*,•,•,*,♦)  is  given  by  (6l).  Formula  (65)  shows 
how  the  joint  distribution  of  in(t,x),  T)(t)  and  Tl(t,x)  is 
relate!  to  the  basic  parameters  of  the  M|G|l  queue. 

The  Limiting  Joint  Distribution.  “Hie  Limiting  joint  distribution 
of  the  th.ee  virtual  waiting  times  is  given  by  its  Laplace- 
Stielejes  transform: 

(66) 

S*(C1’<:2’;3)  =  (1-Xa){1+X  oIC(0»fl+C2+C3,^L+C2,C2+C3»Cl,C3}} 
where  qR^  (-,*,*,*»•,•)  is  given  by  (59). 

Moments  of  the  Limiting  Distribution  of  Basic 
Variables  Uj,  j=0,l,2, 3,4. 

Let  ua  denote: 

(C, 0,0, 0,0) 
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(67)  9  S5(x,8) 

=  sQ  *  wr(*)|R«1)-B'(a3)]+xDL-H(x)3iS,(.2)-Bt-k)] 

(68)  y  s  Y(x,j8)  , 

=  h[X+so-XH(x)h(a3)-x(i-H(x))fitslt)] 

(69)  7  ■  ?(*,s) 

=  h[X-XH(x}fr(Sl)-i(l.H(x))h(s2)] 

‘n*  ^z’i) 

*  V0.*)  ■  a  >  0  , 

vhere  the  functions  n  >  0  are  defined  in  (10).  „ 

follows  from  (62)  that: 


i 


Since  the  calculations  are  lengthy  we  state  the  results  only. 


The  matrices  of  first  and  second  moments  are  respectively 
given  byj 


X  a  a 
x 

# 

X  ax  a 

* 

X  ax  a 


and 


;  b 

I 

I  ax  b  +  Xpxa 

,  2 

i  \  *  2  *  * 

=  :  £  a  b  x  a  a  b  X  a  b  +X0  a 

w  X  XX  XX 

i  2  \Z  *  2  2. 

x  \ b  -  “x b  r  aA b  x  V  *  X3x* 

'  #  2  *  2  #*”  * 

,  x  “x  b  T  °xax  b  V  ax  b  X  “x  “x  b  X  ax  b*’V 

I _ 

'  *  *  . 
where  o^,  ax>  Px  are  defined  inflection  3  and  section  4  and: 

a  *  Uq  *  X0/2(l -X^cr2)  , 

b  -  E  U2  •  4  (l-X3a3)'1[V43X2ae2(l-X2o2)"1]  , 

®  O  A 
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That  is,  a  is  the  steady  state  expected  residual  life  length 
of  the  generation  serving  at  time  t  and  b-a2  its  variance. 

Moments  of  the  limiting  Distribution  of 

3(t»x)>  H(t)  and  Tj(t,x) 

Let  vis  denote: 


3  (t,x)(  !  1  1  0  1  0  ! 


A  s 

|  T'(t)  i  A  *  1  l  1  0  0 

(.  __ 

[  T,(t,x)  '10  10  1 

so  that 

from  (63): 

A  =  A  U 
*** 

(74) 

K  A  -  Me.  3 

(75) 

Ej/l  A')  =  A(E.  H.  o')  A' 

E«  £and  EaUU'  are  given  by  (72)  and  (73),  and  substituting 
these  values  in  (74)  and  (75)  and  simplifying  we  obtain: 

(l+2Xa  )  a, 

(76)  Eo  l\=  (l+>.  t)  & 

(l+2Xa*)  a! 

X 
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(77) 


j  b(l+3Xa  +3X2a2) 
I  ^  x 


+  2X0x  a 


EjM *h  j  b(1+  £  ot+  |  Xax 


!  +  |  X2oax)+X0xa 

j 

I 


l  b(l+  |  Xo+3X2a  o!*) 

i  C  XX 


I 


1 


b(l+Xof+-X2a2) 

+  X  p  a 

b(l+  |  a+  |  Xa*  b(l+3Xa*+3X2ax 
+  |  X2aa*)+X3*  a  +  2  X  3*  a 


The  Limiting  Probability 

p(2(t,x)  <  *n(t)}  as  t  y  « 

fXgtXj)  be  the  joint  distribution  of  Ug  and 
given  t  and  QAi(t,x2,x^)  be  the  probability  that  in  (0,t)  the 
queue  has  never  become  empty  and  that  the  variables  Ug  and 
associated  with  the  time  point  t  satisfy  Ug  <  xg  ,  <  x^, 

given  that  at  t=0  there  were  i  >  1  customers  in  the  queue,  one 
of  who  was  beginning  his  service  at  that  time.  Hien  the 
renewal  argument  as  in  (5)  leads  to: 

pt 

(78)  Ai(t,x2,x3)  =  ^(t.Xg.x^)  +  J  QA1(t-u,x2,x3)dML(u) 

o 

+  PU(t)=0  1  W(o)-i)  U<;x2,x3) 

where : 

U(x2,x3)  =  1  if  x2  >  0  and  x3  >  0  , 

■  0  otherwise 
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Further  let  0(t,t ;;  x^x^)  be  the  joint  distribution  of  Ug  and 
given  t  and  t',  then  as  in  (54,): 


(T9)  0(t,t';x.  ,x  )  =  2  2  e-Xt(l-H(x))-Xt 'h(x) 

j2=°  d3=o 


V 


Similar  to  (57)  we  have: 

a>  as  oo  t 


Jy 


(BO) 


A(t,x2,xJ  -Elf  j  d 0oJ“)(T)dH(v)(t+t,-T) 
’  a  *  *  n=0  v=l  Jo  Jo  0  iv 


(O  (t) 


•  0(t-T,t,;X2,X3) 


V 

v=x  o 


d  0Riv(t)  Fv(t“T’X2’X3) 


where  Q\v(m)  is  defined  in  section  6  and: 

(Bl)  Fv(t-T,x2,x3)  =  |  dH^v^(t+t '-t)  0(t-T,t ';x2,x3) 

o 

(to 

Substituting  (80)  in  (78)  and  applying  Smith's  Key  Renewal 
Theorem  (Theorem  4,  Appendix  D)  we  get  the  limiting  distribution 
A(xg,x3)  of  Ai(t>x2,x3)  as  t-»®,  as  in  (52): 


(02) 

A(x2,x3)=(l-Xa)ju(x2,x3) 


x  A  oRij(+<B)i0Fj(T»x2>x3)dT)  * 


if  1  -  \  a  >  0 


=  0  otherwise 


3  3  C 


Prom  (63)  it  follows  that: 


(83) 


lim  P(T)(t,x)  <  *n(t)}  =  lirn  F(U3  <  Vj 

t  Hr  «*  t  ~>  “ 


l  J  f 


j-1 


(x2)(x3)(t) 


dx2x3Pj^,*xa>x3>,1'r} 


8.  Applications 

Hie  main  objectives  of  a  priority  decision  are  to  reduce 
the  response  time,  to  acknowledge  customer  importance  and  urgency 
of  request  and  to  serve  in  fair  order  and  to  limit  the  length  of 
wait.  Por  the  best  average  performance  the  shortest  service- 
time-next  rule  may  be  just  rig)-  But  under  this  rule  a  steady 
stream  of  shortest  requests  may  delay  a  longest  request  indefinitely. 
The  rul*  proposed  in  this  chapter  is  a  compromise  to  this,  since 
within  each  generation  the  service  request  of  a  customer  with 
long  service  time  is  fulfilled.  Our  model,  of  course,  assumes 
that  the  service  times  of  the  customers  can  be  ordered  before 


hand 
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APPENDIX  A 

A  THEOREM  *N  SUMMATION  OF  SERIES 


Theorem 

For  a  given  positive  integer  k,  the  sum  of  the  infinite 
series : 


(1) 


Z 

n*0 


^(y)v  (yk-D  vy 
>^-1)  yk_1 


for  all  x,  for  all  y/l,  and  for  all  integral  values  of  v  >  0, 
where  [j~]  is  the  greatest  integer  not  exceeding  ^  ,  and 
1  »  are  tne  k-th  roots  of  unity. 

Proof: 

Let  us  denote: 


(2) 

-(x»y)  * 

®  n 

r  x 

E  nr 

n»0  n* 

and 

r80 

(3) 

f(o»y)  ■ 

J 

0 

Then: 

f(*»y)  ■ 

.  t! 

(«0 

I  X- 

n«0  i 

Suppose 

that  0  <  v 

<  k-1, 

n*v1v 

r]> 
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a  k-v-1  ,  i.  Zk-v-1  *  gy  3k- v*l  « 

=  ’pi**  J  j&r**  i  7+r+" 

n°0  ■  n*k-v  i*  n*2k-v  8 


svfe  to  -M1'®''®**  -3 


.k-v  , 
a  -1 


i^'T^rh:  * for  W  •=  '*1  > 

i  (f-i)  »k  (*-i)(s -yk) 


(5) 


_k  v/  k  - »  a 

8  ♦  a  (y  -1)  -  y 

<»-i)  («V) 


Next  we  consider  r  k  <  v  <  (r+l)  k-1,  r*0,l,2,. 

- 

f(«,y)  -  s  ^-^s— 


n*0  a 


•  Z  i 


n-0  8 


n+1 


In  thia  case: 


(6) 


[2£-]k 


,rk  £  *- 

n*0  8 


where  r  *  *  v  -  rk  and  0  <  r  *  <  k-1 
Hence  substitution  of  (5)  in  (6)  yields: 

A  rkr-k  _r k  , x  k-. 

?(.. y)  .  r.  H-*  «r  M .-■»-•■  1-1 

(.•1)  (.*  -  y; 


yrV  ♦  ."rV  -  1)  -  & 

(*  -  1)  (.*  -  y) 


(7) 


,  for  »y|  <  Is', 


Ikk 


To  find  the  inverse  trarsform  we  use  Bateman  (1954) ,  Tables 
of  Integral  Transforms  (p232). 

That  is,  if  f(s)  -  , 

where  P(s)  »  (*-a^)  ...(s-a^),  /  ar  for  i  ^  r  and  Q(s) 

is  a  polynomial  of  degree  <  n-1,  then  the  inverse  transform 
* 

of  f(a)  is  given  by: 


f(*>  -  i,  fJ-j 

o=l  m'  m 


where  P_(a) 
n 


s*«m  * 


Comparing  with  this  we  have  in  (7): 


P(s)  .  («-l)  (.k  .  yk) 


a  (s-1)  (s-u>oy)  (s-o^y)  •••  (s-uj^y) 

so  that  =  1,  o^g  *  u=0,l,. . .  ,k-l 

ai  ^  for  ^  8*nce  yjk  by  assumption.  Where  tao,c^,.  •  • 

are  the  roots  of  ik  *  1  *  0 


«(•>  -  y*[.k  *  .'’"V-D  -  yk] 


Pj^)  ■=  1  -  yk 

W<W>  *  (V1J  >,l"1fv0o!(V“r  ••• 
IVViXWi1  •••  <Wi> 

•  s-  (v1^'1*  -0.1.. ...k-l 
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where  the  simplification  is  obtained  from  the  properties  of 
the  roots  of  the  equation  a*  -  1  =  0. 

QCo^)  -  0 

*  (“^)v  (yk-i)»  ®sOti, . . . ,k-i. 

Hence  using  (8)  the  inverse  transform  f(x,y)  of  f(s,y)  is 
given  by: 


f(x,y) 


k+1  Q(o^) 

m=l  pflJ"°m^ 


x 


1  "I1  %  v 

E  Jo  (V-D  y*'1  * 


This  is  independent  of  r  and  hence  the  result  is  true  for  all 


v  >  0. 
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APPENDIX  B 

PROPERTIES  OF  THE  TABOO  PROBABILITIES 
oQ^(*)  FOR  THE  TANDEM  QUEUE  WITH  ZERO 
SWITCHING 

(The  results  and  proofs  of  this  appendix  closely  follow  Neuts  (l96?>) 

Starting  with  the  seai-Markov  sequence  1^,^,  n  >  0} 
defined  in  (1.3)  we  define  the  taboo  probabilities  oQ^(*) 
as: 

■  8« 

(1)  s  \3M  -  p{5„  ■  i.  *D  <*  !  W  -  ‘5 

and 

0q|^(x)  =  p{T1+...+Tn  <x,  K0=i) 

n  >  1  , 

Let  oq|^(s)  be  the  L.S.T  of  oQ^(x)  and 

(2)  04n^(s>z)  =  z  o0^(8)  zi  » 

Further  ve  denote: 

(3)  °‘u{,) '  it  «4i")(,)  *li0> 


I1*? 


w  o"o3(s)  *  &  o\)(s) 

It  is  seen  that: 

(5)  S  0-y(*)  -  £  [0lin>(a.*)  -  o4n)'8>°)] 

j=l  n=l 

We  define  the  following  sequence  of  functions: 
ao(s,z)  *  z 

(6) 

*n(»>*)  *  \{s,h2[s+X-Xan_1(S,z)]}  ,  n  >  1  , 
where  Y,(*,«)  is  defined  in  (1.12). 

Throughout  this  appendix  we  use  the  following  notations: 

-  h 

<<0’1)  ■  $  C“(M)L 


vj(0,i)  »  ^  \(o,z)l 

3z  z=l 


Lemnft  1 

If  1-X  -  X  a2  >  0  then: 


(7) 


I 

n«i 


and 


]>8 


Proof; 

From  (6)  we  obtain; 

a^(0,l)  =  X  a2  V^(0,1)  VJ0*!)  >  »  >  1  > 

Successive  substitution  yields; 

(9)  an(0»l)  =  [*>  <*2  ao(0>l) 


since  y'(0,1)  =  by  C1*20) 

1  Xoig 

Hence  (7)  follows  for  <  1. 

Similarly  differentiating  (6)  twice  and  simplifying  results 
in  equation  (8). 

Lemma  2 

For  R(s)  >  0  and  i  >  1, 

(1°)  oqi°)(S,z)  =  ^s’z) 

=  a*(s,z)  -  &^_i(8  >&)>  n  >  1, 

Proof. 

From  (2.3)  we  have; 
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(^)  ^(sjz)  =  {Sjhg(  s+X-Xz ) } 

*  a*(s,z) 

Again, 

qi[s,a1(s,z)]  =  Y^Khgts+x-Xa^M)]) 

-  ag(s,z) 
and 

(12)  l  04”’1)(.)  4v[8,h(s,r): 

V— 1 

-  \  Arl)< ■>  *2<  *■*> 

v=l 

-  ,i)]  -  04°'1)(!‘ .0) 

(13)  -  l  ^(S)  «(.,•) 

v=l  v 

=  o^}(3)  a^s»z) 

=  o4n)ts»«i(8»2)3  ■  o4d)(s»0) 

Setting  z  *  0  in  (13)  leads  to: 

(14)  o<4n)(«,0)  =  oq[n'l)C«,a1(s,0)]  -  oqjn_l)(8,0) 
Substitution  of  (12)  and  (l4)  in  (13)  yields: 
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4n+l)(s,z)  =  oq[n'1^[s,a2!(s,2)]  -  0lin”J';[s»a1( 


(n-l)  i 


Successive  similar  substitution  gives: 


(15)  0<4n+l)(s,z)  -  0^i;[B,an(sjZ)]  -  oqi^CB,an.1(s,0)] 

=  V^(s,h2[s+X-Xan(s,z)]} 

-  V^{s,h2[s+x-Xan_1(s,0)]) 

=  an+i(s>z)  '  an(s>°) 

From  (10)  and  (15)  we  have: 

Snh*  *\  =  .V-  1 


,(D 


(1), 


For  n  =  0, 


_q^n;(s,z)  =  a*(s,z)  -  ftn-1(-8»0^  n  -  1> 

oq[°^(s,z)  =  z1 

=  a^(sjz) 


Lemma  3 


Far  s  >  0  and  i  >  15 

(16)  lim  an(s,0)  =  ^m^s)  =  V(s) 

n  ■)  40 

Proof: 

From  lemma  2  we  have: 

.(n) 


;(s,z)  =  an(s,z)  -  Vl^’0^  n  -  l’ 


8,0)] 
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Hence: 

<17)  n-1  ’  aB(s'0) 

The  left  side  of  (17)  is  the  L.S.T.  of  the  probability 
Ajj(x),  where  AK(x)  is  the  probability  that  a  busy  period  with 
one  customer  initially  lasts  for  at  most  N  cycles  of  tasks 
and  has  a  duration  at  most  x.  It  can  be  argued  as  in  Neuts 
(1969): 

<  ak+1(x)  <  1> 

which  implies  that  the  transforms  a^(s,0)  is  increasing  in 
N  for  s  >  0.  Hence  by  Helly-Bray  theorem  (Theoreu  2,  Appendix  D), 
aN(s,°)  converges  to  the  L.S.T.  of  a  probability  mass  function. 
That  is: 

lim  a  (s,0)  =  E  o4n)(B,0) 

N  ®  a  n=l  0  1 

=  o^o(^  =  Y(s)  (1‘35) 

Lemma  4 

If  E  Jn)(s,0)  uP  =  Y(s,o),  |d  <1, 
n=l  0  1 

then: 

E  3  Yi(8>u,)>  1  >  1. 

0=1  0  1 


(18) 
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Proof: 

Analogous  to  the  proof  of  lemma  1.2. 

Lemma  5 

If  R(s)  >  0  and  !<J,  <1,  then  z  -  Y(s,a>)  is  a  root  of 
the  equation: 


(19)  2  =  o)  Y^ts  ,h2(a+X-Xz)]  ,  |z  !  <1, 


Further  z  =  y(s,<b)  is  the  only  root  of  this  equation  in  the  unit 
circle  jz  j  <1  if  R(s)  >  0  and  |o|  <1  or  R(s)  >  0  and  |u)|  <  1 
or  R(s)  >  0,  |u|  <  1  and  1-Xa1-Xa2<0. 

Proof: 

Consider  the  recurrence  relation: 


(20) 

which  gives: 


q(n+l)(.0 

o\o  KS' 


CO 

s  <1 

v=l 


.  (s)  q^(s)  ,  n  >  1, 

lv  0  VO  v  7  — 


os 

E  U> 
n=l 


n+i  0^ri,(.>  -  ®  \  4lv<«>  ’ 

v=l  n=l 


E  d 


i(n}(‘) 

O  VO  v 


=  w  E  4,  (s)  yV(s»w) 

v=l 

(by  lemma  4) 


(21)  =  u>  (q1[s,Y(s,uj)]  -  qlo(s)) 

That  is: 


Y(s,u>) 


CO 

E 

R=1 


n 

UJ 


AnM 


C®* v(  s>yv)3 


(22) 
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=  a)  Y1[s,h2Ce+\-XY(s,uj)]} 

(by  equation  (2.3)) 
which  proves  the  first  part  of  the  lemma. 

The  second  part  follows  from  Rouche's  theorem  (Theorem  3, 
Appendix  D).  For  a  complete  proof  we  refer  to  Takacs  (1962), 
p.  48. 

Lemma  6 

For  R(s)  >0,  2  =  y(b)  is  a  root  of  the  equation: 

(23)  z  =  Y1{s,h2(s+X-Xz)}  ,  jzj  <1  , 

Further  Y(o)  is  the  smallest  positive  real  root  of  the  equation: 

(24)  e  =  y1{o,  h2(x-xe)} 

and  if  1-X  -  X  a2  <  0  then  Q  <  1,  if  1-X  -  X  a2  >  0  then 

0=1. 

Proof: 

Proof  of  the  first  part  is  similar  to  that  of  lemma  5,  by 
taking  uj  =  1. 

For  the  second  part,  the  proof  is  analogous  to  that  of 
lemma  2  in  Neuts  (1969).  For  completeness  we  repeat  that  proof 
here,  since  our  functional  equation  is  different  from  that  in 
Neuts  (1969). 

Consider  the  graphs  of: 

y  =  x  and  y  ■  Y1(0,h2(X-Xx)} 
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and  consider  the  increasing  sequence  of  points  whose  abscissae 
are  an(0,0).  At  the  point  x  =  an(0,0)  we  have: 

y  =  V1{0,hL[X-Xan(0J0)])=  an+1(0,0) 

and  lim  a  (0,0)  a  Y(0)  which  implies  that  Y(o)  is  the  smallest 
n  •  n 

positive  real  root  of  (24). 

If  1-X  i  •  X  a2  <  0  then  from  lemma  5  it  follows  that  9  <  1. 
If  1"X  0|  •  X  a2  >  0,  then  the  graph  of  y  =  Y^OjhgCX-Xx)}  does 
not  intersect  the  line  y  =  x  in  [0,1)  so  that  9  =  1  is  the  only 
root  of  (24). 

Remark : 

From  lemmas  5  and  6  it  follows  that  if  1-X  -  X  a2  >  0 

then: 


(25>  •  v''0)  ■  rr&ffi-;- 

3,  +  2  a,  a,  +  30 
(26)  y'(0)  -  1  1  2  2 


(1-X  0^  -  X  «2) 


T 


(27)  I;  v(0,O)]^0  «  v'(0,l)  -  1  ^  \  a2 

Lemma  7 

If  s  >  0,  0  <  z  <  1  and  l-Xa^-Xag  >  0  then: 


(28) 


4  n\  [l  •  8n(0-0)] 
n=:i  n=i 


X  a„ 

2 

X  a,  -  \  a. 


< 
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Proof: 

The  summands  an(s}z)  -  an(s,0)  is  a  monotonic  increasing 
function  of  z  and  a  decreasing  function  of  a.  Hence  by  setting 
a  ~  0  and  z  =  1  in  the  summands  we  get  the  first  part  of  the 
inequality  (28). 

It  can  be  shown  that  a*(0,z)  >  0  for  all  z  in  [0,1],  so 
that  an(0,z)  is  a  convex  function  in  [0,1]  and  its  graph  lies 
entirely  above  the  tangent  at  z  =  1.  This  tangent  has  an 
intercept : 


\  Ct  ** 

an(0,l>  -  an'(0,l)  -  1  -  (j^) 

where  the  value  of  a#(0,l)  is  taken  from  (9). 

a  a2  n 

Hence  an(0,0)  >  1  -  (■^-~-)  which  proves  the  lemma  completely. 

Limiting  Properties  of  the  Semi -Markov  Sequence.  The  limiting 
properties  of  the  semi-Markov  sequence  defined  in  (1.3)  is 
studied  through  the  following  theorems. 

Theorem  1 

If  lim  P(<>  =j  |  §  =i]  =  P,,  j  >  0,  then  p  =  0  for  all 
n  4  <»  n  0  J  J 

j  if  1-Xa  -  Xo^  <  0.  If  l-Xc^-Xag  >  0  and  0  <  z  <  1,  then 
1  2 


(29) 


B(z)  *  E  H,  zJ  =  1  -  Q  I  [1  -  a  (0,z)] 


j=0 


n=l 


where  pQ  is  given  by: 


t3C) 


i1  =  I  [1  -  &„(O,0)j 


n=0 


and  3^  >  0  for  all  j. 
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Proof: 


The  stationarity  equations  for  the  imbedded  Markov  sequence 
ty  are: 

(31) 


...  ux  irrdGM(u)e-^>M-ma„(v)(T.u) 

J  r-l  r  v=lJo  Ju  ro  J  *  * 


*  ».  .  *  I.  J_  *  £>(u)  ='X(V'U)  3  ‘4w)<m>-  i  ^  o. 


1  tj  <J 

V=1  o  u 


The  first  term  is  obtained  by  considering  r  >  i  customers 
in  unit  1  followed  by  a  1-task  and  a  2-task.  The  second  term 
is  obtained  by  considering  an  idle  period  followed  by  a  1-task 
and  a  2-task.  Equation  (31)  shows  that  all  0.  are  strictly 

tJ 

positive  if  end  only  if  0  is  strictly  positive. 

c 

From  (31)  we  also  obtain: 

BOO  -  2  0r  E  ^l\o)hl(\-\z)  *  0O  E  g^^Ojh^X-Xz) 

r=l  v=l  v=l 


=  E  0r  Y^{o,h2(X-Xz))  +  3o  Y^O.h^X-Xz)} 
r=l 

=  B(Y1(0,h2(x-Xz)})  -  3q[1  -  Y1(0,h2(X-Xz))] 

That  is: 

(32)  B(Y1{0,h2(X-Xz)})  -  B(t)  -  0oCl  -  Y1fO,h2(x-Xz)}] 
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If  l-Xcc^-Xctg  <  0  then  from  lemma  6,  z  =  Y(o)  <  1  is  a  root  of 
the  equation  z  =  Y^[0,h2(X-Xz)}.  Setting  z  =  Y(0)  in 

(32)  we  get  =  0  which  implies  and  implied  by  0^  =  0  for  all 
j.  If  l-Xa^-XOg  >  0,  we  replace  z  in  (32)  by  ar(0,z) 

r  =  0,1,..., n-1  and  add  the  resulting  equations  to  get: 

n 

(33)  B(an(0,z))  -  B(ao(0,z))  =  3Q  I  [1  -  ay(0,z)] 

ts  "  r-1 

Letting  n  00  and  noting  that  an(0,z)-j  1  as  n  -»  ®  for  every 
z  in  [0,1]  we  have: 

(34)  B(l)  -  B(z)  =3  Z  [1  -  an(0,z)] 

n=l 

which  proves  (29)  since  B(l)  =  1.  Finally  equation  (30)  is 
obtained  by  setting  z  =  0  in  (34). 

Theorem  2 

If  0  <  z  <  1  and  l-Xa^-Xcfg  >  0  then, 

(35)  E  -Jfc.^0)  z<]  =  0”1  Z  31 

d-i  0  j-1  3 

Proof: 


From  (5)  we  have: 
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(by  Theorem  2) 


Similarly, 


(41) 


CD  C5 

2  j(j-l)  hl  (0)  =  I  a'(0,l) 
j=l  n=l  n 


16^ 


Substitution  of  lemma  1  in  (4o)  and  (4l)  proves  the  theorem. 


r 
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APPENDIX  C 

To  facilitate  reading  Chapter  IV  we  state  and  prove  the 
following  lemmas  which  ar°  essentially  in  Neute  (1969). 
Starting  with  the  semi-Markov  sequence 

U(Tn),  Tr+1  -  Tn>  n  >  0}  defined  in  (4.1)  we  define  the 
taboo  probabilities  as; 

®  J 


oQif(x)  *  5ij  U(X) 


(1)  o^l)(x)  SQij(x)  =  p^(^n)=0>  Vi-InfxIrtTj-i} 


o*iJ  <x)  -  PfTn  <x,  S(Tn)-d,5(Tv)/0fv-ll...,a-l|?(To).i} 
oqi"'^*)  be  the  t.S.T.  of  0Q?^(x)  and 


,*)  =  L  o^(s)  fz  I  <  1 , 


We  define  the  fallowing  sequence  af  functions: 


ho(s,z)  =  2 


\(s,z)  =  h[«*\  -  \  h^i,,)]  ,  n  >  1 


If2 


Further  we  denote : 


hn(0'1)  '  h  V°'*>] 


Z  =  1 


hn(0'1>  ’  h> 

br.  “7=1 


*>.u .  £  „ 


*(0,,)li 


Lemma  1 

If  1  -  X  a  >  0,  then: 


00 


(5) 


(6) 


£  <(0,D  -  A 


n-0 


£  <(0,1) 

n=0  n 


E  h"'(0,l) 
n=0 


1-Xa 


X23 


(X-Xor)  (l-xV) 


i  r  x 

Ti^TL- 


X3V 


3  x5  a  3 2 


X3a3  (l-X2/)(l-xV) 


;] 


where  a,  3  and  Y  are  defined  in  Chapter  IV. 

The  proof  is  similar  to  that  of  lemma  1  in  Appendix  B. 


Lemma  2 


(7) 


For  R(s)  >  0  and  1  >  1, 
oqi  '*(s>z)  =  <(s,z) 


*>qi  (®>z)  =  hn(8>z)  "  n  >  1 


\ 
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Again  the  proof  of  this  lemma  is  analogous  to  the  proof 
of  lemma  ?,  in  Appendix  B. 

Lemma  3 

If  s  >  0,  0  <  z  <1  and  1  -  \  a  >  0,  then: 


(8) 


£  [h  (s,z)-h  (0,0)] 
n- 1 


<  £ 
n-1 


[i-hn(©,o)] 


< 


1 a 

l-\a 


Frr  the  proof  we  refer  to  the  parallel  proof  of  lemma  7 
in  Arpendix  B. 

F*r  further  properties  of  the  taboo  probabilities  defined 
in  this  appendix  we  refer  to  Neuts  (1969)  where  they  are 
extensively  treated. 


\ 
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APPENDIX  D 

r  SOME  WELL  KNOWN  THEOREMS  USED 

IN  THE  TEXT 

Theorem  1;  ZYGMUND'S  THEOREM 

Let  {Fn(x)}  be  a  sequence  of  distribution  functions  all 
vanishing  for  x  <  0  and  let 
r®  . 

(1)  <?(<«)  =  e1UX  dF  (x),  -  00  <  id  <  00  . 

n  J  n 

o 

-  If  the  functions  0n(ui)  tend  to  a  limit  in  an  interval  around 

u>  =  0,  and  if  the  limiting  function  is  continuous  at  a  =0, 
then  there  i6  a  distribution  function  F(x)  such  that  Fn(;t) 
tends  to  F(x)  at  every  point  of  continuity  of  F(x). 

[Ref:  Zygmund  (1951)] 

Theorem  2:  HELL Y- BRAY  THEOREM 

If  g(x)  is  bounded  and  continuous  when  -  ®  <  x  <  ®  and  the 
sequence  of  distribution  functions  Ffi(x)  converges  to  a 
distribution  function  F(x),  then: 

r  “  r  ^ 

(2)  li®  J  «(x)  d  F  (x)  =  1  g(x)  dF(x) 

ft  CO  _®  _CS 


[Ref:  Loeve  (1963)] 
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Theorem  3:  ROUCHE'S  THEOREM 

If  ffc  )  and  g(z)  are  regular  inside  and  on  a  closed 
contour  C,  and  |g(a)|  <  1  f(z)j  on  C,  then  f(z)  and 
f(z)  +  g(z)  have  the  same  number  of  zeros  inside  C. 

[Ref:  Whittaker  and  Watson  (1952)] 

Theorem  4;  SMITH'S  KEY  RENEWAL  THEOREM 

If  M(t)  is  the  expected  number  of  renewals  in  (0,t], 

Q( • )  is  a  positive  integrable  and  decreasing  function,  then: 


TJ  <30 

J  Q(t-u)  dM(u)  ^  j  Q(u)  du 


where  p,  is  the  mean  renewal  time 
[Ref:  Smith  (1958)] 

Theorem  5:  A  TAU3ERIAN  THEOREM 

If  M(t)  is  non-decreasing  and  such  that 

CD 

p  St 

m(s)  =  |  e”  dM(t)  converges  for  R(s)  >  0,  and  if  for  some 
J  o 

a 

non-negative  number  a,  lim  s  m(s)  =  c,  then: 

6  o 

lim  _1_- 

t  *  t°  i^a+l ) 


[Ref:  Widder  (191*7)] 
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